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ABSTRACT 

The deepest XMM-Newton mosaic map of the central 1.5° of the Galaxy is presented, includ¬ 
ing a total of about 1.5 Ms of EPIC-pn cleaned exposures in the central 15” and about 200 ks 
outside. This compendium presents broad-band X-ray continuum maps, soft X-ray intensity 
maps, a decomposition into spectral components and a comparison of the X-ray maps with 
emission at other wavelengths. Newly-discovered extended features, such as supernova rem¬ 
nants (SNRs), superbubbles and X-ray filaments are reported. We provide an atlas of extended 
features within ±1 degree of Sgr A*. We discover the presence of a coherent X-ray emitting 
region peaking around GO. 1-0.1 and surrounded by the ring of cold, mid-IR-emitting material 
known from previous work as the ’’Radio Arc Bubble” and with the addition of the X-ray data 
now appears to be a candidate superbubble. Sgr As bipolar lobes show sharp edges, suggest¬ 
ing that they could be the remnant, collimated by the circumnuclear disc, of a SN explosion 
that created the recently discovered magnetar, SGR J1745-2900. Soft X-ray features, most 
probably from SNRs, are observed to fill holes in the dust distribution, and to indicate a di¬ 
rect interaction between SN explosions and Galactic center (GC) molecular clouds. We also 
discover warm plasma at high Galactic latitude, showing a sharp edge to its distribution that 
correlates with the location of known radio/mid-IR features such as the ”GC Lobe”. These 
features might be associated with an inhomogeneous hot ’’atmosphere” over the GC, perhaps 
fed by continuous or episodic outflows of mass and energy from the GC region. 

Key words: Galaxy: centre; nucleus; interstellar medium; ISM: supernova remnants; bub¬ 
bles; kinematics and dynamics; X-rays: binaries; diffuse background; ISM; plasmas; methods: 
data analysis; 


1 INTRODUCTION 

At a distance of only ~ 8 kpc, the center of the Milky Way is the 
closest Galactic nucleus, allowing us to directly image, with incom¬ 
parable spatial resolution, the physical processes typical of galactic 
nuclei. The central region of the Galaxy is one of the richest labora¬ 
tories for astrophysics (Genzel et al. 2010; Morris et al. 2012; Ponti 
et al. 2013). Within the inner ~ 200 pc about 3 — 5 x 10^ M© of 
molecular material are concentrated, the so called Central Molec¬ 
ular Zone (CMZ). This corresponds to about 1 % of the molecular 
mass of the entire Galaxy and it is concentrated in a region of about 
~ 10~® of its volume (Morris & Serabyn 1996). In this region 
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many thousands of persistent and transient point-like X-ray sources 
are embedded, such as active stars, bright accreting binary systems 
(and many more quiescent massive bodies) and cataclysmic vari¬ 
ables, which have been beautifully imaged thanks to the superior 
spatial resolution of Chandra (Wang et al. 2002; Muno et al. 2003; 
2009). One of the best jewels in the GC is Sgr A*, the electromag¬ 
netic counterpart of the closest supermassive black hole (BH; Gen¬ 
zel et al. 2010). In addition to this large population of point sources, 
extended X-ray sources, such as supernova remnants, non-thermal 
filaments, pulsar wind nebulae, and massive star clusters, populate 
the GC (Wang et al. 2002). The GC is considered a mini-starburst 
environment, giving us the possibility to study the interaction be¬ 
tween supernova remnants (SNRs) and molecular clouds and the 
impact of massive-and-young star clusters on their surroundings. It 
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allows us to image, in superb detail, the creation and evolution of 
bubbles and superbubbles and the generation of Galactic outflows, 
powered by past starbursts and/or accretion events onto Sgr A*, and 
their impact on the GC environment. 

Warm {kT ~ 1 keV) and hot {kT ~ 6.5 keV) thermal plasma 
emission plus non-thermal hard X-ray emission associated with X- 
ray reflection nebulae (see Ponti et al. 2013 for a review) pervade 
the central region, producing a high background of soft and hard X- 
ray radiation. About 90% (Ebisawa et al. 2001; Wang et al. 2002) 
of the soft X-ray emission appears to be due to a diffuse, patchy 
and thermal component (Bamba et al. 2002) with a temperature 
kT ~ 1 keV, most probably associated with supernova remnants. 
The origin of the hot component is, instead, still highly debated. At 
~ 1.5° from the GC, ~ 80 % of this emission has been resolved 
into point sources (e.g., accreting white dwarfs and coronally ac¬ 
tive stars) by a deep Chandra observation (Revnivtsev et al. 2009). 
Although the intensity of the hot plasma emission increases rapidly 
towards the GC, point sources continue to make a substantial con¬ 
tribution to the observed hard emission (Muno et al. 2004; Heard 
& Warwick 2013a). Additionally, some of the emission may arise 
due to scattering of the radiation from bright X-ray binaries by the 
dense interstellar medium (Sunyaev et al. 1993; Molaro et al. 2014). 
Nevertheless, it is not excluded that a truly diffuse hot-plasma com¬ 
ponent is also present in the GC (Koyama et al. 2009; Uchiyama et 
al. 2013). Such hot plasma would be unbound to the Galaxy and it 
would require a huge energy {E ~ 10^^ erg) and energy loss rate 
of the mass outflow of ~ 10^^ erg s“^, corresponding to a rate of 1 
supernova/yr, to continuously replenish it (Tanaka 2002). However, 
it has recently been proposed that such hot plasma might be trapped 
by the GC magnetic held (Nishiyama et al. 2013). 

Indeed, the magnetic held is thought to be an important in¬ 
gredient of the GC environment. The first high-resolution radio 
images of the Milky Way center (see bottom panel of Fig. 6), re¬ 
vealed the presence of many straight, long (up to ~ 20 — 30 pc) 
and thin (with width <0.1 pc), linearly polarised vertical filaments 
with spectral index consistent with synchrotron radiation (Yusef- 
Zadeh et al. 1984; 1987a,b; Anantharamaiah et al. 1991; Lang et 
al. 1999; LaRosa et al. 2000). These filaments are hypothesized to 
be magnetic flux tubes trapping energetic electrons and therefore 
tracing the diffuse interstellar GC poloidal magnetic held (Morris 
& Yusef-Zadeh 1985; Lang et al. 1999). A staggeringly powerful 
poloidal magnetic held pervading the GC, with a held strength of 
B > 50 fiG, and very possibly B ^ 1 mG, has been inferred 
(Morris 1990; Crocker et al. 2010; Ferriere et al. 2011). The de¬ 
tails of the physical process creating the filaments and energising 
the magnetic held are still debated; however, it appears clear that 
the magnetic filaments are interacting with the ionised surfaces of 
massive molecular clouds. 

Recent far-infrared/sub-millimeter polarization studies of 
thermal dust emission made it possible to probe the direction of 
the interstellar magnetic held inside dense molecular clouds. The 
magnetic held threading GC molecular clouds is measured to be 
parallel to the Galactic plane (Novak et al. 2003; Chuss et al. 2003; 
Nishiyama et al. 2009). Therefore, it appears that the large-scale 
GC magnetic held is poloidal in the diffuse interstellar medium 
and toroidal in dense regions in the plane. If the strength of the dif¬ 
fuse magnetic held is on the high side (B ~ 1 mG) a huge amount 
of magnetic energy, E ~ 10^^ erg, would be stored in the central 
~ 300 pc. This is comparable to the kinetic energy associated with 
the rotation of the gas in the CMZ. Therefore, it is thought to be a 
key player for the GC physics and phenomenology. 

A large scale structure with a possible magnetic origin and 


appearing to be interacting with massive clouds of the CMZ (simi¬ 
lar to the non-thermal filaments) is the Galactic center lobe (GCL). 
The GCL has a limb brightened shell structure in the 10.5 GHz 
map, defined primarily by two spurs (see Fig. 1 and 2 of Law et al. 
2009). The eastern one arises from the location of the GC Radio 
Arc[^ while the second starts from the Sgr C thread. It was pro¬ 
posed that the GCL is produced by channelling of plasma from the 
Galactic plane, induced by energetic GC activity (e.g. episode of 
AGN activity, or a large mass outflow due to the high star forma¬ 
tion rate, etc.; see Law et al. 2011) or from twisting of poloidal 
magnetic held lines by Galactic rotation (Sofue et al. 1984; 1985; 
Uchida 1985; 1990; Shibata 1989). Located at the western limb of 
the GCL is an interesting feature, AFGL 5376 (Uchida et al. 1994), 
an unusually warm, shock heated and extended IR source, thought 
to be associated with the GCL. 

All major X-ray telescopes devoted a significant fraction of 
their time to the study of the GC. Chandra invested several Ms to 
monitor both Sgr A'^’s activity (Baganoff et al. 2001; 2003; Neilsen 
et al. 2013) as well as diffuse soft and hard X-ray emission (Wang 
et al. 2002; Park et al. 2004). Suzaku and Swift also performed large 
observational campaigns to scan the Milky Way center (Koyama et 
al. 2007; Degenaar & Wijnands 2010) and monitor the transients 
in the region (Degenaar et al. 2012). The study of the GC is one 
of the key programs of the NuSTAR mission (Harrison et al. 2013; 
Barriere et al. 2014; Mori et al. 2013). XMM-Newton completed a 
first shallow (~ 30 ks total cleaned exposure in each point) scan 
of the CMZ within a couple of years after launch (see the confer¬ 
ence proceedings: Sakano et al. 2003; 2004; Decourchelle et al. 
2003). A larger amount of time (more than ~ 1.5 Ms) has been 
invested by XMM-Newton on studying the emission properties of 
Sgr A* (Goldwurm et al. 2003; Porquet et al. 2003; 2008; Belanger 
et al. 2005; Trap et al. 2011; Mossoux et al. 2014), focussing on the 
central ~ 15 arcmin, only. Using the XMM-Newton observations 
from the shallow scan of the CMZ together with a number of the 
Sgr A* pointings. Heard & Warwick (2013a,b) have investigated 
the distribution of the X-ray emission within the central region of 
the Galaxy. With the aim of studying the propagation of echoes of 
the past GC activity within the CMZ (Sunyaev et al. 1993; Koyama 
et al. 1996; 2008; Revnivtsev et al. 2004; Muno et al. 2007; Inui 
et al. 2009; Ponti et al. 2010; 2013; Terrier et al. 2010; Nobukawa 
et al. 2011; Capelli et al. 2011; 2012; Clavel et al. 2013; 2014; 
Krivonos et al. 2014), recently, a new deep (with ~ 100 ks ex¬ 
posure at each location) XMM-Newton scan of the CMZ has been 
completed (in fall 2012). We present here the combined images of 
both the new and old XMM-Newton scans, as well as all the XMM- 
Newton observations within the central degree of the Galaxy. 

In ^we present the data reduction process and the key steps to 
produce the GC EPIC mosaic maps. Section 3 introduces the broad 
band X-ray images, discussing the (transient) emission from the 
brightest point sources, the contribution from the foreground emis¬ 
sion, as well as the soft and hard GC diffuse emission. In section 4, 
the narrow band images at the energies of the soft X-ray lines are 
displayed. Section 5 presents a new technique of spectral-imaging 
decomposition of the soft X-ray emission into three physical com¬ 
ponents. Section 6 presents an atlas of all the new and known dif¬ 
fuse features within the surveyed area. Section 7 presents the com¬ 
parison with the distribution of column density of intervening mat- 


^ This is a well known radio feature (see Fig.|^ composed of an array of 
straight, long, thin and linearly polarised vertical filaments, indicating the 
importance of the GC magnetic field. 
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ter. Discussion and conclusions are in § 8 and 9, respectively. Here¬ 
inafter, unless otherwise stated, we will state all locations and po¬ 
sitions in Galactic coordinates and Galactic cardinal points. Errors 
are given at 90 % confidence for one interesting parameter. 


2 DATA REDUCTION AND CLEANING 

The new XMM-Newton CMZ scan has been performed in 2012 
starting on August 30^^ and ending on October 10^^. It comprises 
of 16 XMM-Newton observations all performed with all the EPIC 
instruments in full-frame CCD readout mode with the medium op¬ 
tical blocking filters applied (we refer to Tab. and for more 
details on the instruments set-ups). 

This paper is not limited to the use of the 2012 XMM-Newton 
scan of the CMZ. Instead it is using all XMM-Newton observations 
pointed within 1 degree from Sgr A*. Therefore we combined the 
16 observations of the new CMZ scan with the 14 observations of 
the previous CMZ scan accumulated between 2000 and 2002. We 
also include the 30 observations pointing at Sgr A* and other 49 
observations aimed at studies of different sources in the vicinity of 
Sgr A* (see Tab. an d|9j. 

We performed the analysis of the EPIC data with the version 
13.0.0 of the XMM-Newton Science Analysis System (SAS). Peri¬ 
ods of increased particle background have been removed from the 
data. To perform this, we first selected the Good Time Intervals 
(GTI) starting from the 7-15 keV background light curves, then we 
applied a threshold of 8 and 2.5 cts ks“^ arcmin“^ for EPIC pn and 
EPIC MOS, respectively (see e.g. Haberl et al. 2012). The chosen 
thresholds efficiently cut out all the periods of most extreme activ¬ 
ity of soft proton fiares. We noted, however, that an enhanced, but 
weak, background activity was still present in the data during sev¬ 
eral observations. Because of the non-uniform distribution of the 
GC diffuse emission, lowering the threshold uniformly in all obser¬ 
vations, would result in cutting truly good time intervals in observa¬ 
tions with higher GC diffuse emission. Thus we decided to visually 
inspect the background light curve of all data-sets and select a dif¬ 
ferent threshold for each observation (see Tab. and [^. Such as 
in Haberl et al. (2012), when the data from several detectors were 
available, we combined the GTIs using only common time inter¬ 
vals, otherwise we included GTIs of the single detector. Most of 
the 2012 CMZ scan data were affected by negligible particle fiar- 
ing activity. On the other hand, many of the previous observations 
have been severely affected by soft proton fiares (see the reduction 
in exposure in Tab.[^and[^. 

To prevent infrared, optical and UV photons from bright 
sources in the field of view that would increase the noise and de¬ 
grade the CCD energy scale, the different XMM-Newton observa¬ 
tions have been performed with different filters applied, according 
to the optical-UV brightness of the sources in the field of view (see 
Tab. and [^. In particular, we used the filter wheel closed obser¬ 
vations to remove the internal EPIC background. 


2.1 Images and exposure maps 

Images and exposure maps, corrected for vignetting, have been pro¬ 
duced with an image pixel size of 2” x 2” for each energy band 
(for the definition of all bands, see § |2.2| ). To increase the sky cov¬ 
erage, we selected EPIC-pn events requiring (ELAG & OxfaOOOO) 
= 0, which also includes events in pixels next to bad pixels or 
bad columns. Moreover, we used single to double pixel events. 


EPIC-MOS events were required to have ELAG = 0 and single to 
quadruple-pixel events were allowed. 

Eigure[^ shows the combined EPIC exposure map that covers 
the entire CMZ. Such as done in Sturm et al. (2013), EPIC-MOS 1 
and -MOS2 exposures are weighted by a factor of 0.4 relative to 
EPIC-pn, before being added to the latter, to account for the lower 
effective area. Therefore, the exposure times obtained correspond 
to the equivalent total EPIC-pn exposure time. This allows us to 
obtain a better combination of EPIC-pn and EPIC-MOS data for 
image display purposes. We note, however, that the fluxes can not 
be easily read out directly from these combined images. Therefore, 
the line profiles and the measured fluxes/luminosities are computed 
from the EPIC-pn and each EPIC-MOS map separately and then 
combined (averaged) to obtain a better signal to noise. 

The top panel of Pig. shows that more than 1.5 Ms of clean 
(after cut of time intervals during increased particle background 
activity) exposure time (EPIC-pn equivalent) has been accumulated 
on Sgr A"^ and over ~ 100 — 200 ks are present in each point of the 
CMZ. The few pointings above and below the plane have exposures 
between ~ 15 — 40 ks. Regions with less than 7.2 ks of equivalent 
EPIC-pn exposure have been masked out. 

To check the impact of the bright transients on the images and 
on the physical quantities under investigation, two sets of maps 
have been created. The first series keeps all bright transients and 
point sources, while the second set removes their emission by excis¬ 
ing from the data extended regions including the transients when¬ 
ever they were in outburst (see section[TTJ. The middle and bottom 
panels of Pig. show the exposure maps (computed in the same 
way) for the observations of new and old CMZ scans, separately. 
The maximum exposure times are ~ 190 ks and ~ 45 ks during 
the new and old scan, respectively. 

2.2 Energy bands 

We created images in several energy bands (see Tab. [^. Pigure 
1^ shows the EPIC-pn spectra of the extended emission from sev¬ 
eral regions within the CMZ. In red and black are the spectra from 
the GO. 11-0.11 and Center Superbubble regions, respectively (see 
Pig. Both regions are located within 15 arcmin from Sgr A"^, 
thus they have excellent statistics because of the large exposure. In 
green and blue are the pn spectra of GO. 687-0.146 and Sgr B1 soft, 
respectively, both are located further out, thus having lower expo¬ 
sure. 

We first selected the standard broad energy bands for the con¬ 
tinuum with the softer band being: = 0.5 — 2 keV; the medium 
2 — 4.5 keV and the hard 4.5 — 12 keV (see Tab.[T}. We note that, 
at high energies, the EPIC-pn camera shows strong instrumental 
background emission lines due to Ni Ka (at E ~ 7.47 keV), Cu 
Ka (~ 8.04 keV) and ZnCu (~ 8.63 and 8.87 keV) that strongly 
contribute to the observed X-ray emission in the hard band (see 
Preyberg et al. 2004). To avoid contamination from these strong 
internal background lines, we do not consider (for the EPIC-pn im¬ 
ages) photons in the 1.2-9.2 keV range (see Pig. and Tab.[^. We 
chose these broad energy bands because they are typically used as 
input by the standard XMM-Newton point source detection algo¬ 
rithm and for comparison to other similar surveys of nearby galax¬ 
ies (e.g. M33: Misanovic et al. 2006; Tiillmann et al. 2011; M31: 
Henze et al. 2014; Stiele et al. 2011; EMC: Haberl et al. 1999; 
SMC: Haberl et al. 2012; Sturm et al. 2013). However we note that, 
given the typical GC neutral column density of several 10^^ cm“^, 
the low energy absorption cut-off occurs at the highest energies of 
the standard soft band, making standard broad band RGB images 
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Figure 1. (Top panel) Combined exposure map of all the XMM-Newton EPIC-pn + MOSl + MOS2 observations within one degree from Sgr A*. Such as 
done in Sturm et al. (2013), EPIC-MOSl and -MOS2 exposure is weighted by a factor of 0.4 relative to EPIC-pn to account for the lower effective area. 
The exposure times, thus, correspond to the equivalent total EPIC-pn exposure time. Regions with less than 7.2 ks of equivalent EPIC-pn exposure have been 
masked out. The cleaned EPIC-pn equivalent exposure time is reported in seconds. (Medium panel) Similar exposure map for the observations of the new 
CMZ scan only. (Bottom panel) Similar exposure map for the old CMZ scan only (regions with less than 7.2 ks of equivalent EPIC-pn exposure are included). 
The maximum exposure times are ~ 1.5 Ms, ~ 190 ks and ~ 45 ks during the total, new and old scan, respectively. 


poorly sensitive to column density fluctuations. For this reason we 
define a second set of broad bands, the ”GC continuum bands” (see 
Tab.[^. The first band (E^ = 0.5 — 1.5 keV) contains mainly emis¬ 
sion from foreground sources. The second band (E — 1.5 — 2.6 
keV) is selected in order to contain the low-energy GC neutral ab¬ 
sorption cut-off, thus making it more sensitive to either column 


density or soft gas temperature variations. While the ”GC medium” 
(E; = 2.6 - 4.5 keV) and the ”GC hard” (E; = 4.5 -12 keV) bands 
are similar to the standard broad bands. 

We also selected images at the energies of the soft emission 
lines, such as Si Xlll, S XV, Ar XVII and Ca XIX. To perform contin¬ 
uum subtracted line intensity maps as well as line equivalent width 
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Figure 2. EPIC-pn spectra of the regions: GO. 11-0.11 (red), Center Super¬ 
bubble (black), G0.687-0.146 (green) and SgrBl soft (blue). In dark green 
are the energy bands of the broad GC continuum. In orange (bottom right) 
is shown the part of the hard energy band excluded in order to avoid con¬ 
tamination by Ni Ka, Cu Ka and ZnCu instrumental background emission 
lines. Yellow stripes show the energy bands of the soft and Fe K emission 
lines. Blue stripes indicate the regions selected for the determination of the 
amount of continuum underlying the soft lines. The dotted lines, from top 
to bottom, show the predicted emission of a source with a power-law spec¬ 
trum (with slope r = 1.6) if absorbed by a column density of A^h = 3, 5 
and 9 x 10^^ cm“^, respectively. The blue and orange labels indicate the 
selected broad energy bands for the determination of the continuum under¬ 
lying the Fe K line emission. 


maps, it is essential to measure the level of the continuum under¬ 
lying the emission line. Therefore, we created also several images 
in energy bands on each side of the soft emission lineQ (selecting, 
as far as possible, energy ranges free from line emission; see Fig.[^ 
and Tab.[TJ. 

In the Fe K region we selected two energy bands for the Fe Ka 
and Fe XXV emission. At energies higher than Fe XXV the presence 
of both Fe K/3, Fe XXVI, and of the Fe K edge can give a signifi¬ 
cant contribution. At even higher energies (E ~ 7.5 — 8 keV) the 
contribution from internal background emission line (in the EPIC- 
pn camera) becomes very important, thus we decided to determine 
the continuum emission underlying the Ee K line emission (impor¬ 
tant to determine the Ee K line intensities and equivalent widths) 
through the extrapolation of the continuum red-ward of the Ee K 
lines (see Eig. and Tab. B- The Ee K line emission and its vari¬ 
ations will be the prime scientific focus of two future publication 
(Ponti et al. in prep.; Soldi et al. in prep.; see also Ponti et al. 2014; 
Soldi et al. 2014) and will not be discussed here any further. 

All images were exposure corrected and, to remove readout 
streaks, the images from EPIC-pn were corrected for out-of-time 
events. Noisy CCDs in the MOS data (Kuntz & Snowden 2008) 
have been searched with the SAS task emtaglenoise and removed 
from the mosaic images. 


^ The energy band red-ward of the Si xiii line extends only down to 
1.65 keV, because of the presence of the strong Al Ka background emission 
line at E ~ 1.49 keV (Freyberg et al. 2004). 


Standard continuum bands: 


Soft 

0.5-2 

Medium 

2-4.5 

Hardf 

4.5-12 



GC continuum bands: 


Fore 

GC Soft 

GC Medium GC Hardf 


0.5-1.5 

1.5-2.6 

2.6-4.5 4.5-12 



Soft emission lines: 


Si XIII 

S XV 

Ar XVII Ca xix 


1.80-1.93 

2.35-2.56 

3.03-3.22 3.78-3.99 


Continuum subtraction soft emission lines: 


Red-Si 

Si-S 

S-Ar Ar-Ca 

Blue-Ca 

1.65-1.77 

2.1-2.3 

2.7-2.97 3.27-3.73 

4.07-4.5 



Fe K lines: 


Fe Ka 

Fe XXV 



6.3-6.5 

6.62-6.8 




Continuum subtraction Fe K: 


CFeK 

CsFeK 

CmFeK ChFeK 



soft 

medium hard 


5-6.1 

4.0-4.7 

4.7-5.4 5.4-6.1 



Table 1. Energy bands used for each of the different continuum, and nar¬ 
row line images. Also shown are the energy bands used to determine the 
continuum underlying the line emission. Units are in keV. Several energy 
bands, at lower energies compared to the FeK lines, have been computed 
to determine the best continuum subtraction for the FeK lines. fTo avoid 
contribution from the strong internal detector background emission lines, 
present in the EPIC-pn camera (such as: Ni Ka, Cu Ka and ZnCu), we do 
not consider photons in the 7.2-9.2 keV from this instrument (on the other 
hand, we consider such photons detected in the EPIC-MOS cameras). 


2.2.1 Stray-light rejection 

Because the GC region is crowded with many bright (transient) X- 
ray sources, several observations, including the new XMM-Newton 
scan, are badly affected by stray-light (see Fig. [^. Stray-light is 
produced by photons from sources located outside of the XMM- 
Newton EPIC instrument’s fields of view and singly refiected by 
the mirror hyperbolas, thus creating concentric arc-like structures 
in the detector plane (see XMM-Newton user handbook). The stray- 
light contribution is small (the effective collecting area for stray- 
light is less than ~ 0.2 % of the effective on-axis area), but a very 
bright source can have an important impact up to ~ 1.4 deg outside 
the field of view. 

Analogously to the removal of bright transients, we masked 
the strongest stray-light artefacts in the images of individual obser¬ 
vations. In most cases, affected regions are covered by other un¬ 
affected observations, thus leaving no features in the final mosaic 
map. To remove a stray-light artefact, we defined a rough region 
including the artefact and an individual cut-off value of the surface 
brightness. Using this cut-off, we created a mask from an image of 
this region in the total energy band that has been smoothed with a 
Gaussian kernel with a FWHM of 10^^ beforehand. This mask was 
multiplied by all images and exposure maps of this observation. 

2.2.2 Adaptive smoothing 

All images have been smoothed separately using the SAS tool AS- 
MOOTH. To prevent different smoothing patterns from introducing 
colour artefacts in RGB images, we adaptively smoothed all en¬ 
ergy bands in such images with the same smoothing template. For 
the broadband XMM-Newton continuum RGB images (see Fig. 
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Figure 3. Standard broad energy band (see Tab. [T]and Fig. i RGB mosaic image of all XMM-Newton observations within one degree of Sgr A* (see Tab.f 
This represents the deepest X-ray view of the CMZ region with exposure higher than 0.2 Ms along the Galactic disc and 1.5 Ms in the center (see Fig.[ 
X-ray emission from X-ray binaries, star clusters, supernova remnants, bubbles and superbubbles, HII regions, PWNs, non-thermal filaments, nearby X-ray 
active stars, the supermassive BH Sgr A* and many other features are observed (see Fig. [^and|^. The detector background has been subtracted and adaptive 
smoothing applied. Residual features and holes generated by correction of the stray light from GX 3+1 are visible (see also Fig. 1) at Galactic latitudes between 
I 1.2° and / ~ 1.4° and latitudes h —0.2° and b ~ 0.4°. © ©00 RAS, MNRAS 000,^?? 
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Figure 4. Standard broad energy band RGB image of the XMM-Newton CMZ scan performed in 2012. The CMZ is observed with a uniform exposure (see 
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Figure 5. Finding chart. The brightest X-ray point sources (all X-ray binaries) are labelled in white (see Tab.|^. In red the positions of some star clusters 
are reported, which are placed either in the GC or along the spiral arms of the Galaxy (see Tab.|^. With yellow dashed lines the location of some molecular 
complexes are shown. See www.mpe.mpg.de/heg/gc/ for a higher resolution version of this hgure. 


g0|||7l[^and[nj, we required a minimum signal-to-noise ratio 
of 6 in the 0.5-12 keV energy band image (i.e., the sum of the three 
energy bands composing the RGB image), as well as the standard 
minimum and maximum size of the smoothing Gaussian kernel of 
10" (full width half maximum) and 200", respectively. The signal 
to noise ratio at each pixel is defined as the value at that pixel di¬ 
vided by its standard deviation and the adaptive smoothing that we 
applied is making the signal to noise ratio as close at possible to 
6, therefore fainter or less exposed areas are more smoothed than 
brighter or better exposed regions. For the narrower-hand soft-line 
images (see Fig. and [^, we use the S xv map as a 

template, requiring a minimum signal-to-noise ratio of 4 and the 
same standard minimum and maximum of the smoothing kernel. 
The same smoothing kernel is then applied to all the other bands of 
the RGB images. 


2.2.3 Internal particle background subtraction 

Unless otherwise specified, internal particle background has been 
removed from each broad-band image. Following Haberl et al. 
(2012) we first create, for each selected energy band, both the total 
emission and \hQ filter wheel closed images. We then re-normalise 
and subtract iht filter wheel closed images from the total emission 
images. The filter wheel closed image re-normalisation factor is 
computed by equating, for each instrument, the number of pho¬ 
tons in the unexposed corners of the detector to that in the filter 
wheel closed images (see Haberl et al. 2012 for more details). This 
procedure is reliable and accurate for reasonably long exposures 
{t ~ 5 — 10 ks). For this reason datasets with total clean EPIC-pn 


exposure shorter than 5 ks have not been considered in this analy¬ 
sis. 

2.2.4 Continuum subtraction 

To subtract the continuum emission from an emission-line image, 
we define a narrow band (B) containing the line, typically sand¬ 
wiched by two nearby but generally wider energy bands (A and 
C) that are dominated by continuum emission. Under the assump¬ 
tions that the emission in the A and C bands is dominated by the 
continuum and that the continuum emission can be described by a 
simple power-law, we could in principle determine the intensity of 
the continuum for each pixel of the band B image. Indeed, using 
the fiuxes in A and C bands, it is possible to derive the continuum 
parameters (spectral index F and intensity). However, this requires 
the solution of non-linear equations. Therefore, we prefer to im¬ 
plement a different technique based on interpolation. Using Xspec 
we simulate, for power-law spectral indices going from V = 0.3 
to 3.6, the ratio between the observed fiux (e.g., number of pho¬ 
tons measured) in the continuum in bands A and C, compared to 
the simulated continuum fiux in band B (e.g. Nb/(2 x Na) and 
Nb/(2 X Nc)). We record these ratios and then plot them as a 
function of the hardness ratios {Nc — Na)/{Nc + Na)), which 
is a proxy for the spectral index V. We then find the best-fitting 
linear relationship between these values, thus determining Couab 
and LiriAB that are then allowing us to measure the continuum 
emission underlying the line emission in band B {Nb) from the in¬ 
tensity in band A {Na) and the hardness ratio {Nb = 2 x Na x 
[Couab + LiuAB x {Nc — Na) /{Nc + Na)] - To reduce the un¬ 
certainties, we perform the same procedure for band C, determining 
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Figure 6. Finding charts. (Top panel) Broadband X-ray continuum image. White ellipses show the position and size of known, radio-detected supernova 
remnants. Cyan ellipses indicate the position and size of bright diffuse X-ray emission possibly associated with supernova remnants that lack a clear radio 
counterpart (or such in the case of G359.12-0.05 that show X-ray emission significantly displaced from the radio emission associated to the radio remnant 
G359.07-0.02). The magenta ellipses show the location and dimension of some bright HII regions, while the red ellipses indicate some of the largest non- 
thermal filaments detected in radio (see Tab.[^. Blue ellipses show some PWN and the yellow dashed ellipses show the regions used to accumulate the spectra 
shown in Fig.[^ (Bottom panel) 90-cm radio image of the CMZ region obtained with the VLA (courtesy of LaRosa et al. 2000). For display purposes the radio 
supernova remnants are shown with green ellipses. The other regions have the same colour code as the top panel. The image shows the radio flux (Jy beam“^ 
unit). See www.mpe.mpg.de/heg/gc/ for a higher resolution version of these figures. 
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Galactic longitude 

Figure 7. Finding chart. Zoom of the central ~ 10 arcmin of the Milky Way as seen by XMM-Newton (same energy bands and color scheme as in Fig.[^. The 
position of Sgr A* is indicated by the blue cross. The red ellipses show the position and extent of filamentary and diffuse X-ray emission features associated 
with, e.g., non-thermal filaments (Tab.[^. The magenta dashed ellipses show the location and extension of the 20 pc bipolar X-ray lobes. The black dashed 
ellipses indicate the position and orientation of the CND. See WWW.mpe.mpg.de/heg/GC/ for a higher resolution version of this figure. 


ConcB and LincB- We then average these values obtaining, for 
each pixel: Nb = NaX [Couab + LiriAB x {Nc - Na)/{Nc + 
Na)] + Nc X [ConcB + LincB x [Nc — Na)/{Nc + ^a)]. 
We finally subtract this continuum emission image from the total 
emission image B to determine the line intensity map. 


3 THE XMM-Newton BROADBAND VIEW OF THE 
GALACTIC CENTRE 

Figure shows the broad energy band mosaic image of all exist¬ 
ing XMM-Newton observations within 1 degree from Sgr A*. Fig- 
ure|^shows the Galactic centre image obtained only with data from 
the 2012 XMM-Newton campaign. At the GC distance of 7.8 kpc 
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Figure 8. Finding chart. Chandra RGB image of all the ACIS-I observations pointed at Sgr A* (see Clavel et al. 2013 for data reduction and details on 
the image production). The red, green and blue images show the GC soft (1.5-2.6 keV), GC medium (2.6-4.5 keV) and GC hard (4.5-8 keV) energy bands, 
respectively. The same regions displayed in Fig.^are evidenced here. See www.mpe.mpg.de/heg/gc/ for a higher resolution version of this figure. 


(Boehle et al. 2015), 1 arcmin corresponds to 2.3 pc, 10 pc subtends 
~ 4.3’ and ~ 0.2° corresponds to 28 pc. In red, green and blue, the 
soft (0.5-2 keV), medium (2-4.5 keV) and hard (4.5-12 keV) con¬ 
tinuum bands are shown, respectively. Hundreds of point sources 
and strong diffuse emission are clearly observed in the map. These 
point sources are characterised by a wide variety of colours, rang¬ 
ing from distinctive red to dark blue. 


3.1 Bright and transient GC point sources during the new 
(2012) XMM-Newton CMZ scan 

Many X-ray point sources are clearly visible in Fig. |^and|^ A de¬ 
tailed catalogue of the properties of all the detected point sources is 
beyond the scope of this paper. Here we briefly describe the bright¬ 
est GC sources detected by XMM-Newton and the X-ray transients 
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Figure 9. Zoom of the finding chart displayed in Fig.|^ 


in the field of the 2012 scan (see Degenaar et al. 2012 for a com¬ 
pendium of previously noted transients). 

The brightest GC point source of the 2012 XMM-Newton 
CMZ scan is IE 1743.1-2843 (Porquet et al. 2003; Del Santo et 
al. 2006) a persistently accreting neutron star binary detected at an 
observed fiux level of F 2 -iokeV ~ 1-1 x 10“^° ergcm“^ s“^ (im¬ 
plying an unabsorbed flux of F 2 -iokeV,unab ~ 2.6 x 10“^° erg 
cm“^ s“^; Nh ~ 2 x 10^^ cm“^; obsid: 0694641201). During the 
2012 XMM-Newton campaign we also detected an outburst from a 
new, very faint X-ray transient that we name XMMU J174505.3- 
291445. The source has a typical quiescent luminosity at or be¬ 
low Lx ~ 10^^ erg s“^, but on 2012 August 31^^ (during obsid 
0694640201) it was observed to go into outburst for about ~ 2 hr 
and to reach a peak X-ray luminosity of Lx ~ 10^^ erg s“^. The 
detailed spectral and multiwavelength analysis of this new transient 
will be presented in a separate paper (Soldi et al. in prep.; but see 
also Soldi et al. 2014). 

Another faint X-ray transient, XMM J174457-2850.3, is de¬ 
tected in two 2012 observations; obsid: 0694641101 - 0694640301. 
The observed 2-10 keV flux is 1.1 zb 0.3 and 2.9 zb 0.6 x 
10“^^ erg cm“^ s“^, respectively. A power-law fit to the spectrum 
with the photon index fixed to the value reported by Sakano et al. 
(2005) yields a column density of Nh = (1.4 zb 0.4) x 10^^ cm“^. 
This source was discovered in 2001 by Sakano et al. (2005) who re¬ 


ported a tentative detection of an X-ray pulsation of ~ 5 s, during 
the ~ 25 ks XMM-Newton observation. Both a visual inspection 
and timing analysis of the X-ray light curve show no evidence for 
bursts and/or dips. However, even considering the 4 times longer 
exposure of the new data, we cannot exclude or confirm the ~ 5 s 
periodic modulation because of the lower flux observed. In fact, 
during the 2001 XMM-Newton observation (obsid: 0112972101) 
XMM J174457-2850.3 had a flux about 10 — 40 times higher 
(~ 45 X 10“^^ erg cm“^ s“^) than in 2012 (in quiescence XMM 
J174457-2850.3 has atypical 2-10 keV fiux lower than 0.2 x 10“^^ 
erg cm“^ s“^). 

Only upper limits are measured for the other well-known 
X-ray transients within the field of view. The two bursters 
GRS 1741.9-2853 and AX J1745.6-2901 (Sakano et al. 2002; 
Trap et al. 2009; Ponti et al. 2014; 2015) have fiux limits of 
F2-i0keV < 2 X 10“^^ erg cm“^ s“^ and F2-iokeV < 10“^^ 
erg cm“^ s“^ (obsid: 0694641101, 0694640301), respectively. 
Closer to Sgr A*, we find an upper limit on the 2-10 keV fiux 
of F 2 -iokeV < 5 X 10“^^ erg cm“^ s“^ toward three other 
sources: CXOGC J174540.0-290031, the low-mass X-ray binary 
showing X-ray eclipses (Porquet et al. 2005; Muno et al. 2005), 
CXOGC J174540.0-290005 (Koch et al. 2014), and the magne- 
tar discovered on April 25, 2013 (Degenaar et al. 2013; Dwelly 
& Ponti 2013; Mori et al. 2013; Rea et al. 2013; Kaspi et al. 2014; 
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Coti-Zelati et al. 2015), located at distances from Sgr A* of only 
~ 2.9, ~ 23 and ~ 2.4 arcsec, respectively. 

Finally we observe that both XMMU J174554.4-285456, the 
faint transient with a possible pulsation period of about 172 s (Por- 
quet et al. 2005), and SAX J1747.7-2853, the bursting (showing 
also superbursts) X-ray transient (Wijnands et al. 2002; Natalucci 
et al. 2004; Werner et al. 2004), have flux limits of F 2 -iokeV < 
2 X 10“^^ erg cm“^ s“\ 

Bright sources outside the 2012 scan 

Three very bright sources are outside the held of view during 
the 2012 CMZ scan, however they imprint their presence through 
bright stray-light arcs. The arc features between and south of the 
Sgr A and C complexes (/ - 359.6 - 359.9°, b - -0.15 - 0.4°) 
testify that the bright X-ray burster lA 1742-294 (Belanger et al. 
2006; Kuulkers et al. 2007) was active during the 2012 XMM- 
Newton campaign. The very bright arcs east of the Sgr D complex 
(obsid: 0694641601) are most probably produced by the very bright 
neutron star low-mass X-ray binary GX 3-1-1 (Piraino et al. 2012) lo¬ 
cated about 1.18° northeast of the arc^ On the far west edge of the 
2012 scan a brightening is observed. This is due to IE 1740.7-2942 
(Castro et al. 2013; Reynolds & Miller 2010; Natalucci et al. 2014), 
a bright and persistent accreting microquasar, at only ~ 1.5 arcmin 
from the edge of the 2012 held of view (see Fig. [^. The lack of 
straylight south of the Sgr B region suggests that the BH candidate 
IGR J17497-2821 (Soldi et al. 2006; Paizis et al. 2009) was in qui¬ 
escence during these observations. Two bright X-ray bursters have 
been active during the 2003 XMM-Newton observation pointed to 
the pulsar wind nebula called The Mouse, i.e., SLX 1744-299 and 
SLX 1744-300 (Mori et al. 2005). 

3.2 Very soft emission: Foreground emission 

Despite the presence of distinctively soft (red) point sources, Fig.[^ 
shows no strong, diffuse, very soft X-ray emission. This is mainly 
due to the very high column density of neutral hydrogen toward the 
GC (with typical values in the range Nh ~ 3 — 9 x 10^^ cm“^; 
see also Almost no Galactic centre X-ray radiation reaches us 
below E^l.3, 1.7 or 2.3 keV for column density values ofNn — 
3, 5, or 9 X 10^^ cm“^, respectively (see Fig. |^. The majority 
of the ’’red” sources present in the 0.5-1.5 keV band are point-like 
and are associated with foreground active stars characterised by an 
unabsorbed soft X-ray spectrum. 

Two clearly extended and soft X-ray emitting sources are 
present in Fig.[^ These correspond to Sh2-10 and Sh2-17 (Wang et 
al. 2002; Dutra et al. 2003; Law et al. 2004; Fukuoka et al. 2009), 
two stellar clusters located in one of the Milky Way spiral arms and 
thus characterised by a lower column density of absorbing material, 
consequently appearing stronger in the 1-2.5 keV range (visible in 
Fig.with orange colours). 

3.3 Soft and hard GC emission 

Galactic centre radiation with energies above ~ 2 — 3 keV can 
typically reach us and be detected (in green and blue in Fig.[^. GC 
sources with a signiflcant continuum component (e.g. power-law or 

^ This region is covered only by the observations of the 2012 XMM-Newton 
scan, therefore the removal of the stray-light arcs generates regions with null 
exposures in the final mosaic maps (e.g., Fig. [^and|^. 


Bremstrahlung), such as observed from most GC point sources, the 
GC stellar clusters (e.g. the Arches, the Quintuplet and the Cen¬ 
tral cluster) as well as some supernova remnants (such as SNR 
GO.9-1-0.1 and Sgr A East) appear with a bright light blue colour. 

Colour gradients confirm the presence of at least two compo¬ 
nents of the diffuse emission, each having a different spatial distri¬ 
bution (see Eigs. [^and|^. One component dominates the emission 
in the soft and medium energy bands, thus appearing with a dis¬ 
tinctively green colour. Its distribution appears to be very patchy, 
peaking typically at the position of known supernovae remnants. 
Another, harder component appears with a dark blue colour (Eig. 
1^. This harder emission is known to consist of at least two separate 
contributions. One, which is associated with intense high-ionisation 
Ee K lines, is smoothly distributed and peaks right at the GC; it is 
likely produced by faint point sources (Muno et al. 2004; Revnivt- 
sev et al. 2009; Heard & Warwick 2013a). The other, which is as¬ 
sociated with neutral Ee K emission lines, has a patchy distribution 
peaking at the position of molecular cloud complexes; it is likely 
due to an ensemble of X-ray reflection nebulae (see Ponti et al. 
2013 for a review). We also note that the Galactic plane emission 
is dominated by dark blue colours (Eig. |^, while regions located 
at 5>0.2° and — 0.35° have a significantly greener colour. We 
address this in more detail in §|5]0and|^ 


4 SOFT LINE EMISSION 

Eigurej^shows the spectra of the diffuse emission from the regions 
marked in magenta in Fig.|^ The ~ 1.5 to ~ 5 keV band shows 
strong, narrow emission lines, the strongest of which are Si XIII, S 
XV, Ar XVII and Ca XIX. This line emission, as well as the underly¬ 
ing continuum and the intra-line emission, are typically well fitted 
by a thermal model (e.g. APEC in XSPEC) with temperatures in the 
range 0.6 — 1.5 keV (Kaneda et al. 1997; Tanaka et al. 2000; Muno 
et al. 2004; Nobukawa et al. 2010; Heard & Warwick 2013b). At 
higher energies a power-law component with intense Ee XXV and 
Ee XXVI lines is also observed over the entire GC region. Addi¬ 
tionally, neutral Ee Ka and K/3 lines are also observed. The neutral 
Ee K emission lines are associated with different processes, there¬ 
fore they will be the focus of separate publications. 


4.1 RGB images of soft emission lines 

The top and bottom panels of Eig. show the line (continuum 
non-subtracted) RGB image and the inter-line continuum RGB im¬ 
age (see caption of Eig. and Tab. for more details). We note 
that the soft X-ray line image shows very strong colour gradients 
(less dramatic colour variations are observed in the continuum im¬ 
age). In particular, the sources DSl (the core of Sgr D), the west¬ 
ern part of Sgr B1 (i.e., GO.52-0.046, GO.570-0.001), Sgr C, as 
well as the Chimney above it, all have a distinctively green-blue 
colour, while G359.12-0.05, G359.10-0.5, G359.79-0.26, G359.73- 
0.35 and the entire G359.77-0.09 superbubble are characterised by 
orange-brown colours. GO. 1-0.1, the Radio Arc, the arched fila¬ 
ments (see Eig. of Lang et al. 2002), GO.224-0.032, and GO.40-0.02 
are also characterised by red-brown colours, however here a grada¬ 
tion of white and green is also present (please refer to Tab. and 
l^and Eigs. and for the positions of the regions listed here). 
The lobes of Sgr A appear with a whiter colour than the surround¬ 
ings. In addition, we observe bright red-brown emission along two 
broad, linear ridges having relatively sharp edges to the northwest 
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Figure 10. (Top panel) Soft lines (continuum unsubtracted) RGB image of the CMZ. The Si xiii line emission is shown in red, S xv in green and Ar xvii+ Ca 
XIX in blue. (Bottom panel) RGB image of the energy bands between soft emission lines. The Si-S band emission is shown in red (see Tab.[T]for a definition of 
the energy bands), S-Ar in green and Ar-Ca + Blue-Ca in blue. The diffuse emission in this inter-line continuum is very similar to the soft emission line one, 
suggesting that the same process is producing both the lines and the majority of the soft X-ray continuum. The colour variations within the map are modulated 
primarily by abundance variations (top), temperature of the emitting plasma, continuum shape and absorption. 


and northeast of Sgr A*. This latter feature is discussed in detail in 
section [8^ 

In spite of the fact that these images show different compo¬ 
nents (one being dominated by emission lines, the other by contin¬ 
uum emission), they are remarkably similar. No clear diffuse emis¬ 
sion component is present in one and absent from the other image. 
This indicates that most of the diffuse soft X-ray continuum and 
line emission are, indeed, produced by the same process. In ad¬ 
dition, the differences in the ratio between photons emitted in the 
lines and in the continuum can add valuable information for under¬ 
standing the radiative mechanism. In fact, such differences could be 
due, for example, to different cosmic abundances and/or variations 


in the relative contributions of various thermal and nonthermal radi¬ 
ation mechanisms. In order to better highlight these differences, we 
map the sum of the interline continua in the same image (see cap¬ 
tion of Fig. [TTJand Tab.[^. As expected, none of the point sources is 
a strong soft line emitter (they in fact appear brighter in the interline 
image). We also note that the intense soft X-ray emitting regions in 
the Galactic plane, such as Sgr D, Sgr Bl, Sgr C, the Chimney and 
G359.9-0.125 are characterised by distinctively orange-red colours, 
indicating they are strong line emitters. 

For an alternative perspective, the bottom panel of Fig. m 
shows the Si Xlll + Si-S bands in red, S-Ar + S XV + Ar-Ca in 
green, and CFeK in blue. These energy bands are chosen to high- 
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Figure 11. (Top panel) RGB image composed of summed, continuum-subtracted line emission (Si xiii + S xv + Ar xvii + Ca xix) in red, the sum of the 
interline continua (Si-S + S-Ar + Ar-Ca + Blue-Ca) in green and the CFeK emission in blue (see Tab.[T}. (Bottom panel) RGB image, in red the Si xiii + 
Si-S emission, in green the S-Ar + S xv + Ar-Ca and in blue the CFeK emission. 


light any energy dependence in the soft X-ray emission that could 
be due to column density variations of the obscuring matter or tem¬ 
perature fluctuations of the emitting gas. In fact, the softer energy 
bands (Si Xlll + Si-S) will be more affected by absorption or low 
temperature plasma emission compared to the medium (S-Ar + S 
XV + Ar-Ca) or high energy bands. We defer the detailed discus¬ 
sion of the features present in these images to the discussion of the 
various physical components presented in §[^and subsections. 

4.2 Continuum subtracted soft emission line maps and 
profiles 

From top to bottom, the panels of Fig. show the continuum- 
subtracted Si XIII, S XV, Ar XVII, Ca XIX intensity maps. Although 


the continuum subtraction procedure should naturally remove the 
emission from the line-free point sources, small fluctuations in the 
continuum subtraction, in the case of the brightest sources, some¬ 
times leave signiflcant residuals. For this reason, we have masked 
out the brightest point sources in our computation of these maps. 


The different curves of Fig.[^show the continuum-subtracted 
emission profiles (integrated over latitude from the magenta rectan¬ 
gle in Fig.|12|) for the individual soft emission lines. The same four 
line profiles are compared in Fig. with similar profiles in which 
the contribution of specific bright structures has been removed. 
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Figure 12. From top to bottom, continuum subtracted Si xiii, S xv, Ar xvii, Ca xix intensity maps of all the stacked AMM-iVewron observations of the CMZ. 


5 SPECTRAL DECOMPOSITION 

In order to better trace the relative contributions of the diffuse 
thermal (soft and hot) and non-thermal components, we have per¬ 
formed a simple component separation using a list of images de¬ 
picting various energy bands. We use a total of 17 energy bands: 
11 for the continuunjj and 6 for the lines (tracing Si XIII, S XV, 
Ar XVII, Ca XIX, Fe Ka and Fe XXV, see Tab.[T]). This treatment of 
the data allows us to be more confident about the spectral decom¬ 
position, e.g. compared to single RGB maps, retaining most of the 
morphological information on sufficiently large scales (i.e. beyond 
few arcmin scales). 

The general assumption is that the emission at any position 


The eleven continuum energy bands used are: 1.0-1.5 keV; 1.5-1.8 keV; 
2.0-2.35 keV; 2.55-3.05 keV; 3.25-3.75 keV; 3.95-4.70 keV; 4.70-5.40 
keV; 5.40-6.30 keV; 6.50-6.60 keV; 6.80-7.80 keV and 8.20-9.50 keV. 


can be represented by the linear sum of three main components, 
namely i) a soft plasma with a temperature of 1 keV (Kaneda et al. 
1997; Bamba et al. 2002); ii) a hot plasma of temperature 6.5 keV; 
and iii) a non-thermal component modeled by an absorbed power- 
law plus a neutral, narrow iron line (with 1 keV equivalent width), 
that are subject to an additional absorbing column of Nh = 10^^ 
cm“^. All three components are also absorbed by gas in front of the 
GC region and both thermal plasmas are modelled using the APE 
model in XSpec. The resulting model is therefore PHABS (APEC + 
APEC -f PHABS (POWERLAW -I- Gauss))) and has only three free 
parameters: the relative normalizations of the three components. 
The hot plasma component represents the emission associated with 
faint unresolved point sources, whose cumulative spectrum is well 
described by a thermal spectrum (Revnivtsev et al. 2009) plus a 
possibly truly diffuse hot plasma component (Koyama et al. 2007). 
The spectral index of the non-thermal component is assumed to be 
r = 2, consistent with the values measured through the combined 
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Figure 13. Longitudinal intensity profiles of the Si xiii (red), S xv (green), 

Ar XVII (blue) and Ca xix (violet) emission lines, integrated over Galactic 
latitude within the magenta rectangular region shown in Fig.[T^ 

spectral fits of XMM-Newton spectra with higher energy data (e.g., 
INTEGRAL and/or NuSTAR\ Terrier et al. 2010; Mori et al. 2015; 
Zhang et al. 2015). 

The strongest assumptions in this approach are that the emis¬ 
sion can be represented everywhere with these three components. 
This obviously fails on bright point sources or on regions where 
the emission is much hotter (e.g. Sgr East). For the soft compo¬ 
nents the even stronger assumption is that absorption to the GC is 
assumed to be uniform over the CMZ at a value of N/f = 6 x 10^^ 
cm^^ (Sakano et al. 2002; Ryu et al. 2009), putting aside absorption 
in the GC region itself. Clear column density modulations are ob¬ 
served towards different lines of sight (e.g. Ryu et al. 2009; Ryu et 
al. 2013). We tested significantly different column densities (up to 
Nh = 1.5 X 10^^ cm“^ characteristic of several GC sources; see 
e.g. Baganoff et al. 2003; Rea et al. 2013; Ponti et al. 2015). We 
found that if the soft plasma normalization is significantly modi¬ 
fied, the overall morphology is consistent. We tested various values 
of the other parameters (spectral index or temperatures) and did not 
find strong effects on the soft plasma morphology or normalization. 

We first produced counts, exposure and background maps 
for each observation and each instrument. Background was ob¬ 
tained from cal-closed datasets distributed in the ESA^calibration 
database. For each observation and instrument, an average RMF is 
computed as well as the un-vignetted ARF. For each instrument, 
individual observation images were reprojected using the final im¬ 
age astrometry and then combined to compose a mosaic. Average 
ARF and RMF for each instrument were obtained with the FTOOL^ 
routines ADDRMF and ADDARF. 

For each pixel of the final maps, we fit the measured num¬ 
bers of counts in all the energy bands and instruments with a model 
consisting of the three aforementioned components as well as the 
background events number and the Out-of-Time (OoT) events for 
the EPIC-pn camera. The free parameters are the normalization of 
each individual component. We apply Cash statistics (Cash 1979) 
to take into account the low statistics in each pixel. This analysis 

^ http://xmm 2 .esac.esa.int/external/xmm_sw_cal/background/epic_esas.shtml 
® https://heasarc.gsfc.nasa.gov/ftools/ftools_menu.html 


Nh 

( 1022 ) 

F 2 - 4 . 5 /FO 

F 4 . 5 - 10 /FO 

0.01 

1 

1 

1 

0.776 

0.968 

3 

0.488 

0.912 

5 

0.322 

0.857 

7 

0.223 

0.809 

10 

0.137 

0.737 

15 

0.0698 

0.636 

30 

0.0156 

0.414 

50 

0.0033 

0.244 

70 

0.0008 

0.149 

100 

1.2 X 10-4 

0.0754 

150 

6.2 X 10-^ 

0.0273 


Table 5. Expected ratio of the obscured flux to the un-obscured flux for 
different values of the column density of obscuring material. A thermally 
emitting gas with temperature of kT = 1 keV (phabs*apec model) is as¬ 
sumed in the computations. The predicted observed flux in both the 2 — 4.5 
and 4.5 — 10 keV bands (F 2 - 4.5 and ^ 4 . 5 - 10 ) is computed and compared 
to the respective un-absorbed (FO) flux. 


allows us to perform a rough spectral decomposition, better sep¬ 
arating the spectral emission components, although retaining the 
maximum spatial resolution. 

Figure [^presents the map of the normalisation (in units of 
10“^ times the APEC normalisation) of the soft thermal emission 
component. The normalisation of the soft thermal component has a 
distribution similar to the one traced by the soft lines and the con¬ 
tinuum (Fig.[^[^[TT]and[^. Enhanced high-latitude soft plasma 
emission is observed. The white dashed lines show the position of 
two sharp edges in the distribution of this high latitude emission 
(see also Fig.|^[^and[^. The white solid line shows the edge of 
the region having more than 7.2 ks of exposure (see Fig.[^. 


6 AN ATLAS OF DIFFUSE X-RAY EMITTING 
FEATURES 

The patchy and non-uniform distribution of the diffuse emission 
makes the recognition of the shape, the border and connection of 
the different structures and components difficult. Occasionally, dif¬ 
ferent works report the same X-ray feature with different names 
and shapes and, in extreme cases, the same X-ray emitting feature 
is associated with different larger scale complexes. 

In Tab. and we report all the new X-ray features 
discussed in this paper, plus many GC features presented in 
previous works. The main purpose of these tables is to pro¬ 
vide a first step towards the building of an atlas of dif¬ 
fuse X-ray emitting GC features. The table is available online 
at: www.MPE.MPG. de/HEG/GC/AtlasGCdiffuseX-ray and 
will be updated, should the authors be notified of missing extended 
features. This exercise is clearly prone to incompleteness and defi¬ 
ciencies, however we believe this might help in providing a clearer 
and more systematic picture of the diffuse X-ray emission from 
the GC region. The spatial location and size of all these features is 
shown in the finding charts in Fig. |^and|^ 
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Figure 14. Map of the normalisation of the soft thermal gas component (in units of 10“"^ times the APEC normalisation). The white lines indicate the extent of 
the survey having more than 7.2 ks exposure. The white dashed lines show the position of the two sharp edges in the distribution of the high latitude plasma. 
Some bright point sources (i.e., IE 1743.1-2843, AX J1745.6-2901, IE 1740.7-2942, GRS 1741.9-2853) have been removed, thereby producing artificial holes 
in the maps at their respective locations. 


7 THE FOREGROUND COLUMN DENSITY 

Given the high column densities of neutral or weakly ionized ma¬ 
terial absorbing the soft X-ray radiation, it is important to estimate 
the effects of X-ray obscuration. For example, a molecular complex 
having a column density of Nh ~ 10^^ cm“^, such as the Sgr B2 
core, would completely obscure the radiation below about 4 keV, if 
placed in front of the GC; see Fig.|^ 

To calculate the effects of absorption of the X-ray emission, 
we computed the flux generated by a thermally emitting plasma 
with temperature of kT = 1 keV (using a PHABS*APEC model), in 
both the 2 — 4.5 and 4.5 — 10 keV bands, after being absorbed by a 
given column density of neutral material (see also Fig.[^. For each 
column density explored, we report in Tab. the ratio of the ob¬ 
served flux (F 2 - 4.5 and F 4 . 5 - 10 ) over the respective un-absorbed 
(FO) flux. We note that the hard X-ray band starts to be affected 
(corresponding to flux reductions up to a factor of 2) for column 
densities up to Nh ~ 3 x 10^^ cm“^ while it is heavily affected 
(flux reduction of a factor of 10 or more) for Nh ~ 10^^ cm“^ 
or higher (see Tab.[^. At lower energies, the obscuration effect is 
even more pronounced. Already, for Nh ~ 3 x 10^^ cm“^, the 
observed flux is less than half and for Nh ^ 5 x 10 “ cm"^ it is 
less than 0.1 % of its un-obscured flux. This indicates that the softer 
band is expected to be heavily affected by absorption. 

7.1 Column density distribution 

The top panel of Fig. shows the neutral Hydrogen column den¬ 
sity distribution as derived from dust emission (Molinari et al. 
20110 The image shows the Nh distribution in logarithmic scale 
in the range Nh = 4.5 x 10^^ — 3.8 x 10^^ cm“^. This total column 

^ We do not show the entire CMZ, because of the limited coverage of the 
Herschel dust emission map (Molinari et al. 2011). 


density estimated from the dust has large uncertainties that can be 
mainly ascribed to the uncertainty associated with the dnst-to-Nn 
ratio. In particular, the column densities shown in this map appear 
to be systematically larger than what is measured with other meth¬ 
ods. For example, the column density of the GO. 11-0.11 massive 
cloud is estimated to be Nh ~ 5 — 6 x 10^^ cm“^ in this map, 
while Amo-Baladron et al. (2009) measure Nh ~ 2 x 10^^ cm“^, 
through a detailed modelling of the molecular line emission. The 
core of Sgr B2 is estimated by Molinari et al. (2011) to have 
Nh ~ 3 X 10^^ cm“^, while modelling of the X-ray emission 
(Terrier et al. 2010) suggests Nh ~ 7 x 10^^ cm“^, more than an 
order of magnitude lower. Moreover, the average column densities 
of GO.40-0.02, GO.52-0.046, GO.57-0.018 and a fourth region (the 
magenta ellipse in Fig. are estimated to be Nh ~ 4 x 10^^, 
4 X 10^^, 1.2 X 10^^ and 1.5 x 10^^ cm“^, respectively, from 
the dust map, while they are measured to be in the range Nh ~ 
7 — 10 X 10^^ cm“^, from modelling of the X-ray emission. There¬ 
fore, the total normalisation of the Nh map built from the dust 
emission appears to be overestimated. However, the method em¬ 
ployed to produce it does not suffer from self-absorption, so it is 
presumably giving unbiased relative Nh ratios. 


7.2 X-ray emission modulated by absorption 

The bottom panel of Fig.[^shows the X-ray map with the column- 
density contours overlaid for comparison. We observe that, as ex¬ 
pected, no soft X-ray emission is observed toward the central part 
of the most massive molecular cores. In particular we observe de¬ 
pressed X-ray emission from: i) the Sgr B2 nucleus and its enve¬ 
lope (with Nh > 7 X 10^^ cm“^); ii) the almost perfect coin¬ 
cidence between the hole in soft X-ray emission east of G0.224- 
0.032 (see Fig. and and the shape of the so-called ’’Brick” 
molecular cloud, M0.25-i-0.01 (see Fig.[^ Clark et al. 2013); hi) 
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Galactic lonaitude 



Galajclic lonqitude 


Figure 15. Top panel: Neutral Hydrogen column density distribution as derived from dust emission (Molinari et al. 2011). The image shows the Nh distribu¬ 
tion in logarithmic scale from Nh = 4.5 x 10^^ up to 3.8 x 10^^ cm“^. The green, magenta and white contour levels correspond to Nh = 1.5 x 10^^, 
7 X 10^^ and 1.5 x 10^"^ cm“^, respectively. Bottom panel: X-ray continuum RGB map (Fig.[^ with the column-density contours overlaid. 
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BRIGHT AND TRANSIENT POINT SOURCES 

Source name 

Coordinates^; 

Fluxf 

References 

within the 2012 CMZ scan 

IE 1743.1-2843 

0.2608,-0.0287 

110 

90,92,93 

XMMU J174505.3-291445 

359.6756,-0.0634 

14 

94 

XMMU J174457-2850.3 

0.0076,-0.1743 

0.3 

90,95 

GRS 1741.9-2853 

359.9528,+0.1202 

< 0.02 

90,59,97 

AX J1745.6-2901 

359.9203,-0.0420 

< 0.1 

59,90,91,124 

CXOGC J174540.0-290031 

359.9435,-0.0465 

< 5 

98,99 

SGR J1745-2900 

359.9441,-0.0468 

< 5 

101,102,103,104 

XMMU J174554.4-285456 

0.0506,-0.0429 

< 0.2 

98 

SAX J1747.7-2853 

0.2073,-0.2385 

< 0.2 

90,105,106,107 

CXOCG J174540.0-290005 

359.9497,-0.04269 

< 5 

100 

within the total GC scan 

IE 1740.7-2942 

359.1160,-0.1057 


111,112,113 

lA 1742-294 

359.5590,-0.3882 


90,108,109 

IGRJ17497-2821 

0.9532,-0.4528 


90,114,115 

GRO J1744-28 

0.0445,+0.3015 


90 

XMMU J174654.1-291542 

359.8675,-0.4086 


90 

XMMU J174554.4-285456 

359.1268,-0.3143 


84,85,86 

SLX 1744-299 

359.2961,-0.8892 


37,59,87,88,89 

SLX 1744-300 

359.2565,-0.9111 


37,59,87,88,89 


Table 2. List of bright and transient point sources during the 2012 XMM-Newton scan as well as bright point sources observed in all scans of the region (see 
Fig.[^. To avoid exessive crowding around Sgr A*, CXOCG J174540.0-290005 and SGR J1745-2900 are not shown. fFluxes are given in units of 10“^^ erg 
cm“^ s“^ and correspond to the mean flux observed during the 2012 XMM-Newton scan of the CMZ. ^Coordinates are in Galactic format. References: (1) 
Wang et al. 2006a; (2) Yusef-Zadeh et al. 2002; (3) Capelli et al. 2011; (4) Tatischeff et al. 2012; (5) Sakano et al. 2003; (6) Habibi et al. 2013; (7) Habibi et 
al. 2014; (8) Krivonos et al. 2014; (9) Clavel et al. 2014; (10) Dutra et al. 2003; (11) Law et al. 2004; (12) Fukuoka et al. 2009; (13) Wang et al. 2002a; (14) 
Tsuru et al. 2009; (15) Mori et al. 2008; (16) Mori et al. 2009; (17) Heard & Warwick 2013a; (18) Maeda et al. 2002; (19) Park et al. 2005; (20) Koyama et al. 
2007a; (21) Kassim & Frail 1996; (22) Nobukawa et al. 2008; (23) Senda et al. 2002; (24) Renaud et al. 2006; (25) Mereghetti et al. 1998; (26) Gaensler et al. 
2001; (27) Porquet et al. 2003a; (28) Aharonian et al. 2005; (29) Dubner et al. 2008; (30) Nobukawa et al. 2009; (31) Sawada et al. 2009; (32) Morris et al. 
2003; (33) Morris et al. 2004; (34) Markoff et al. 2010; (35) Zhang et al. 2014; (36) Nynka et al. 2013; (37) Gaensler et al. 2004; (38) Pedlar et al. 1989; (39) 
Cotera et al. 1996; (40) Figer et al. 1999; (41) Johnson et al. 2009; (42) Lu et al. 2008; (43) Lu et al. 2003; (44) Yusef-Zadeh et al. 2005; (45) Baganoff et al. 
2003; (46) Ho et al. 1985; (47) Bamba et al. 2002; (48) LaRosa et al. 2000; (49) Morris & Yusef-Zadeh 1985; (50) Lang et al. 1999; (51) Anantharamaiah et 
al. 1991; (52) Yusef-Zadeh & Morris 1987a; (53) Yusef-Zadeh & Morris 1987b; (54) Yusef-Zadeh & Morris 1987c; (55) Muno et al. 2008; (56) Uchida et al. 
1992; (57) Predehl & Kulkarni 1995; (58) Senda et al. 2003; (59) Sakano et al. 2002; (60) Coil et al. 2000; (61) Murakami 2002; (62) Yusef-Zadeh et al. 2007; 
(63) Dutra & Bica 2000; (64) Zoglauer et al. 2014; (65) Koyama et al. 2007b; (66) Nakashima et al. 2010; (67) Downes & Maxwell 1966; (68) Tanaka et al. 
2009; (69) Tanaka et al. 2007; (70) Wang et al. 2006b; (71) Wang et al. 2002b; (72) Phillips & Marquez-Lugo 2010; (73) Hewitt et al. 2008; (74) Reich & 
Fuerst 1984; (75) Gray 1994; (76) Roy & Bhatnagar 2006; (77) Marquez-Lugo & Phillips 2010; (78) Borkowski et al. 2013; (79) Yamauchi et al. 2014; (80) 
Inui et al. 2009; (81) Green 2014; (82) Yusef-Zadeh et al. 2004; (83) Nord et al. 2004; (84) Uchiyama et al. 2011; (85) Heinke et al. 2009; (86) Muno et al. 
2006; (87) Mori et al. 2005; (88) Skinner et al. 1990; (89) Pavlinski et al. 1994; (90) Degenaar et al. 2012; (91) Ponti et al. 2014; (92) Porquet et al. 2003b; 
(93) Del Santo et al. 2006; (94) Soldi et al. 2014; (95) Sakano et al. 2005; (97) Trap et al. 2009; (98) Porquet et al. 2005a; (99) Muno et al. 2005b; (100) Kock 
et al. 2014; (101) Degenaar et al. 2013; (102) Dwelly & Ponti 2013; (103) Rea et al. 2013; (104) Kaspi et al. 2014; (105) Wijnands et al. 2002; (106) Natalucci 
et al. 2004; (107) Werner et al. 2004; (108) Belanger et al. 2006; (109) Kuulkers et al. 2007; (110) Piraino et al. 2012; (111) Castro et al. 2013; (112) Reynolds 
& Miller 2010; (113) Natalucci et al. 2014; (114) Soldi et al. 2006; (115) Paizis et al. 2009 (116) Lu et al. 2013; (117) Do et al. 2013; (118) Yelda et al. 2014; 
(119) Bamba et al. 2000; (120) Bamba et al. 2009; (121) Ohnishi et al. 2011; (122) Zhao et al. 2013; (123) Hales et al. 2009; (124) Ponti et al. 2015. 


the core of the Sgr C comple}0 iv) the regions around DB-58 and 
at Galactic position / ~ 0.2,6 ~ —0.48° also appear to have darker 
colours and, once again, it is possible to find molecular complexes 
(MO.Ol 8-1-0.126 and MO.20-0.48) covering roughly the same region 
(see Fig.j^an d[Tg. All these clouds are characterised by very high 
column densities N/f >3 — 7 x 10^^ cm“^ and they most proba¬ 
bly lie in front of Sgr A"^ and of most of the GC (e.g., according to 
the twisted ring model of Molinari et al. 2011). Therefore, they are 
absorbing the GC’s extended, soft X-ray emission. 

All this evidence suggests that at least the most massive clouds 

® At this location a sharp transition in the soft X-ray emission, with an 
arc-like shape, is observed. This is spatially coincident to the edge of a very 
dense core of dust, suggesting that the modulation in the soft X-ray emission 
is induced by obscuration by the molecular cloud. 


located in front of the GC do actually modulate (obscure) the soft 
X-ray emission. However, fluctuations in column densities cannot 
be the only cause for the observed distribution of soft X-ray emis¬ 
sion for two reasons. First, we do observe only weak soft X-ray 
emission along some lines of sight having a low column density 
of molecular material (such as around the Sgr C and Sgr D com¬ 
plexes and south of the Sgr B1 region). Second, we do detect in¬ 
tense (among the brightest) soft X-ray emission from several re¬ 
gions such as the Sgr A complex and the cores of the Sgr C and 
Sgr D complexes, where some of the highest column-density clouds 
are found. In particular, in Sgr A very intense soft X-ray emission is 
observed along the line of sight toward the 50 km s“^, the Bridge 
(Ponti et al. 2010) and the GO. 11-0.11 clouds, some of the high¬ 
est column density clouds in the CMZ. Although this might be 
explained by placing these clouds on the far side of the CMZ, it 
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ATLAS OF DIFFUSE X-RAY EMITTING FEATURES 

Name 

Other name or 

Coordinates 

Size 

References 


associated features 

(1, b) 

arcmin 


STAR CLUSTERS: 

Central cluster 


359.9442, -0.046 

0.33 

45,116,117,118 

Quintuplet 


0.1604, -0.0591 

0.5 

1,63,11 

Arches 

G0.12-t0.02 

0.1217, 0.0188 

0.7 

1,2,3,4,5,6,7,8,9,39,40,11 

Sh2-10bt 

DB-6 

0.3072,-0.2000 

1.92 

10,11,12,63,11 

Sh2-17bt 

DB-58 

0.0013, 0.1588 

1.65 

13,63,11 

DB-05bt 

GO.33-0.18 

0.31 -0.19 

0.4 

22,63,11 

SNR and Super-bubbles candidates: 

G359.0-0.9ttt G358.5-0.9 - G359.1-0.9 

359.03,-0.96 

26 X 20 

X-R 48,51,75,76,81,119,120 

G359.07-0.02 

G359.0-0.0 

359.07,-0.02 

22 X 10 

R 14,48,51,66 


G359.12-0.05 

359.12,-0.05 

24 X 16 

X66 

G359.10-0.5ttt 


359.10,-0.51 

22 X 22 

X-R 37,48,51,56,74,75,81,120,121 

G359.41-0.12 


359.41,-0.12 

3.5 X 5.0 

X 14 

Chimney bt 


359.46,-t0.04 

6.8 X 2.3 

X 14 

G359.73-0.35tt 


359.73,-0.35 

4 

X58 

G359.77-0.09 

Superbubble 

359.77,-0.09 

20 X 16 

X 15,16,17,58 


G359.9-0.125 

359.84,-0.14 

15 X 3 

X 15,16,17,58 


G359.79-0.26^ 

359.79,-0.26 

8 X 5.2 

X 15,16,17,58 

20 pc lobes 

G0.0-0.16t^ 

0.00,-0.16 

359.94, -0.04 

5.88 

X This work 

X-R 32,33,34,17 

G359.92-0.09t 

Parachute - G359.93-0.07 

359.93,-0.09 

1 

R 35,38,43,47,58,60,61 

Sgr A East 

GO.O-tO.O 

359.963, -0.053 

3.2 X 2.5 

X-R 5,18,19,20,48,75,81 

GO. 1-0.1 

Arc Bubble 

0.109,-0.108 

13.6 X 11 

X This work 


G0.13,-0.11b 

0.13,-0.12 

3X3 

X 17 

GO.224-0.032 


0.224,-0.032 

2.3 X 4.6 

X This work 

G0.30-t0.04 

G0.3-t0.0 

0.34,-t0.045 

14 X 8.8 

R 21,48,51,81,82 

GO.42-0.04 

G0.34-t0.05 

G0.33-t0.04 

Suzaku J1746.4-2835.4 

0.40,-0.02 

4.7 X 7.4 

X22 

GO.52-0.046 

GO.40-0.02 

0.519,-0.0460 

2.4 X 5.1 

This work 

GO.57-0.001 


0.57,-0.001 

1.5 X 2.9 

This work 

G0.57-0.018t 

CXO 1174702.6-282733 

0.570,-0.018 

0.2 

X 23,24,58,59,68,80 

G0.61-t0.01t 

Suzaku J1747.0-2824.5 

0.61,-tO.Ol 

2.2 X 4.8 

X 22,65,79 

G0.9-t01^ 

SNR0.9-t0.1 

0.867,-t0.073 

7.6 X 7.2 

R 25,26,27,28,29,48,75,81,82 

DSl 

Gl.2-0.0 

1.17,-tO.OO 

3.4 X 6.9 

X31 

Sgr D SNR 

G 1.02-0.18 

1.02,-0.17 

10 X 8.0 

R 30,31,48,51,75,77,81,82 

Gl.4-0.1 

Gl.05-0.15 

Gl.05-0.1 

Gl.0-0.1 

1.4,-0.10 

10 X 10 

R 73,81,82 


Table 3. Atlas of diffuse X-ray emitting features. The first two columns in the table indicate the name primarily used in this work to refer to the feature as well 
as the other names used in the previous literature. The third and fourth columns show the coordinates of each feature as well as its approximate projected size. 
Finally, the fifth column provides references to selected works discussing the feature. For convenience, we report in Tab.[^all the references ordered according 
to the numbering used in this table. The other names column shows the different designations used in previous literature. In the case of bubbles, these features 
are not necessarily referring to the same structure but to features forming the bubble candidate. fPossibly due to a thermal hlament. JThe interpretation as a 
SNR is probably obsolete, ft Most probably a foreground feature. ^ This feature appears to be part of the superbubble G359.77-0.09. New extended X-ray 
feature, possibly part of the superbubble G359.77-0.09. b This feature appears to be part of the Arc bubble, o Possibly connected to GO.61+0.01. ^ X-ray 
emission primarily non-thermal, therefore it appears also in the next table, ft New extended X-ray feature, possibly part of the superbubble G359.77-0.09. bf 
The low X-ray absorption towards these star clusters indicate that they are located in front of the GC region, fft Because of the low X-ray absorption column 
density (Nh ~ 2 x 10^^ cm“^) this is most probably a foreground source (Bamba et al. 2000; 2009). bf The Chimney is most probably either part of a large 
scale structure (see §8.7) or an outflow from G359.41-0.12 (Tsuru et al. 2009), therefore most probably it is not a separate SNR. 


appears that these regions are characterised by truly enhanced soft 
X-ray emission (see e.g. §8.4, 8.5 and 8.6). We defer the detailed 
disentangling of these effects to an elaborate spectral study of these 
regions. 


8 DISCUSSION 

In the process of systematically analysing all XMM-Newton ob¬ 
servations of the central degrees of the Galaxy, we have discov¬ 
ered several new extended features and have produced an atlas of 
known, extended soft X-ray features. Here we discuss their general 
properties and investigate the origin/existence of several specific 
features. 


8.1 General properties 


To compute the total observed (absorbed) fiux from the CMZ 
we first mask out the emission from the brightest binaries, by 
excluding a circle with a radius of: 1' around SAX J1747.7- 
2853; 1.5' for AX J1745.6-2901 and GRS1741.9-2853; 2' for 
XMMU J174445.5-295044; 2.5' for IE 1743.1-2843; 3.5' for 
IGR J17497-2821 and IE 1740.7-2942 (see white circles in Fig. 
1^. We then compute the total observed count rate from two boxes, 
one with a size of 1.5° x 0.35° (I x b) and centered on Sgr A* and 
one with a bigger size of 2.08° x 0.413° centered on I = 0.232°, 
b = 0.080°. We then measure the total count rate within these re- 
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ATLAS OF DIFFUSE X-RAY EMITTING FEATURES 

Name 

Other name or 

Coordinates 

Size 

References 


associate features 

(1, b) 

arcsec 


Radio and X-ray filaments and PWN candidates: 




Snake 

G359.15-0.2 

359.15,-0.17 

312 X 54 

R48 

G539.40-0.08 


359.40,-0.08 

27.5 X 5.1 

X41 

G359.43-0.14 


359.43,-0.14 

21.4 X 3.9 

X41 

Sgr C Thread 


359.45,-0.01 

500 X 42 

R 48,51 

Ripple filament 

G359.54-r0.18 

359.548,-r0.177 

320 X 55 

R 43,44,48,51 

G359.55+0.16 

X-ray thread 

Suzaku J174400-2913 

359.55,-r0.16 

56.1 X 8.0 

X 13,41,42,43,79 

Crescent 

G359.79-r0.17 

Curved filament 

359.791,-r0.16 

300 X 74 

R 63,50,51 

Pelican 

G359.85-r0.47 

359.859,-r0.426 

300 X 54 

R 48,50,51 

Cane 

G359.87-r0.44 

G359.85-r0.39 

359.87,-r0.44 

420 X 50 

R48 

Sgr A-E 

G359.889-0.081-wisp 

XMM J174540-2904.5 

G359.89-0.08 

359.889,-0.081 

20 X 5 

RX 5,35,41,42,43,44,50,55 

G359.897-0.023 


359.897,-0.023 

6.4 X 4 

X55 

G359.899-0.065 

Sgr A-F 

G359.90-0.06 

359.899,-0.065 

6.5 X 2.5 

X 42,44,55 

G359.904-0.047 


359.904,-0.047 

6.5 X 3 

X55 

G359.915-0.061 


359.915,-0.061 

7X2 

X55 

G359.91-1.03 


359.919,-1.033 

138 X 36 

R48 

G359.921-0.030 

F7 

359.921,-0.030 

7.5 X 3 

X 42,55 

G359.921-0.052 


359.921,-0.052 

5.5 X 2 

X55 

The Mouse 

G359.23-0.82 

359.30-0.82 

156 X 108 

PWN F 37,48,57,123 

G359.925-0.051 


359.925,-0.051 

8 X 2.2 

X55 

G359.933-0.037 

F2 

359.934,-0.0372 

12 X 3 

X 41,42,55 

G359.933-0.039 

FI 

359.933,-0.039 

5X2 

X 42,55 

G359.941-0.029 


359.941,-0.029 

6X2 

X 41,55 Stellar wind 

G359.942-0.045 


359.942,-0.045 

5X3 

X55 

G359.944-0.052 


359.944,-0.052 

9 X 1.5 

X 41,55 

G359.945-0.044 


359.945,-0.044 

6 X 2.5 

X 41,1,42,55 PWN 

G359.95-0.04 


359.950,-0.043 

10 X 4 

X 55,70 PWN 

G359.956-0.052 


359.956,-0.052 

4 X 2.5 

X55 

G359.959-0.027 

F5 

359.959,-0.027 

9X3 

X 41,42,55 

Southern thread 

G359.96-r0.09 

359.96-r0.09 

359.96,-rO.ll 

500 X 40 

R 48,50,51 

G359.962-0.062 


359.962,-0.062 

5.5 X 3.5 

X55 

G359.964-0.053 

F3 

359.964,-0.053 

16 X 3.5 

X 41,42,45,55 PWN 

G359.965-0.056 

F4 

359.965,-0.056 

9X3 

X 42,55 

G359.969-0.033 


359.969,-0.033 

5X2 

X55 

G359.970-0.009 

F8 

359.970,-0.009 

10 X 2.5 

X 41,42,55 PWN 

G359.971-0.038 

F6 

359.971,-0.038 

16 X 8 

X 41,42,55 PWN 

G359.974-0.000 

F9 

359.974,-0.000 

7X2 

X42 

G359.977-0.076 


359.977,-0.076 

6X4 

X55 

Cannonball 

J174545.5-285829 

359.983,-0.0459 

30 X 15 

X-R PWN 36,122 

G359.983-0.040 


359.983,-0.040 

6.5 X 4.5 

X 42,55 

G359.98-0.il 


359.979,-0.110 

Streak 

R50 

GO.007-0.014 

GO.008-0.015 

0.008,-0.015 

11 X 3.5 

X 41,55 

GO.014-0.054 


0.014,-0.054 

18 X 14 

X55 

GO.017-0.044 

MC2 

0.017,-0.044 

15 X 4 

X 41,42 FeKa 

G0.02-r0.04 


0.0219,-r0.044 

Streak 

R50 

G0.021-0.051 


0.021,-0.051 

15 X 12 

X55 

GO.029-0.08 


0.029,-0.08 

29 X 18 

X55 

GO.032-0.056 

GO.029-0.06 - FIO 

GO.03-0.06 

0.0324,-0.0554 

35 X 6 

FeKa 41,42,55 PWN 

GO.039-0.077 


0.039,-0.077 

22 X 15 

X55 

G0.062-r0.010 

G0.06-r0.06 

0.062,-r0.010 

40 X 25 

R 50,55 

Northern thread 

G0.09-r0.17 

G0.08-r0.15 

0.09,-r0.17 

714 X 48 

R 48,49,50,51 

GO.097-0.131 


0.097,-0.131 

70 X 50 

X55 

Radio Arc 

GCRA 

GO. 16-0.15 

0.167,-0.07 

1690 X 145 

R 38,48,49,50,51 

GO.l 16-0.111 


0.116,-0.111 

50 X 40 

X55 

GO. 13-0.11 


0.13,-0.11 

55 X 12 

71,17,41,42 PWN 

GO. 15-0.07 

Steep spectrum of Radio Arc 

0.138,-0.077 


R50 

XMM JO. 173-0.413 

GO. 17-0.42 

0.173,-0.413 

180 X 18 

X This work 


S5 

0.17,-0.42 

912R82 


GO.223-0.012 


0.223,-0.012 

50 X 20 

X41 

G0.57-0.018t 

CXO J174702.6-282733 

0.57,-0.0180 

0.33 

X 23,24,58,59,68,79,80 

G0.61-r0.01t 


0.61,-rO.Ol 

132 X 288 

X 22,65,79 

G0.9-r0.1 


0.9,-rO.l 


PWN 


Table 4. Atlas of diffuse X-ray emitting features. This table has the same structure as Tab.[^ Because of possible misplacements between the peak emission 
of the radio and X-ray counterparts of filaments, SNR, PWN and other diffuse structures (generally related to the different ages of the population of electrons 
traced at radio and X-ray bands), when available we give the best X-ray position (following the preference: Chandra, XMM-Newton, Suzaku), otherwise we 
state the radio position. We cite the literature results separating thes between the X-ray (X) from the radio (R) detections. fPossibly part of a young supernova 
remnant. For convenience, we report in Tab. I^all the references ordered according to the numbering used in this table. 
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Figure 16. XMM-Newton image in the 2-12 keV band showing the new 
X-ray filament located south of the Radio Arc. The filament is located at 
I ~ 0.173° and b ~ —0.413° and appears as a thin (< 0.15 arcmin) and 
long (~ 2.7 arcmin) filament running along the north-south direction, as 
the Radio Arc. The brightest source in the image is SAX J1747.7-2853. 

gions and convert it into a flu}|^assuming that the soft X-ray diffuse 
emission is dominated by a thermally emitting plasma with a tem¬ 
perature of /cT = 1 keV (Kaneda et al. 1997; Bamba et al. 2002). 
In particular, we assumed an APEC emission component with tem¬ 
perature kT =1.08 keV, absorbed by a column density of neutral 
material ofNn = 6 x 10^^ and with Solar abundance. 

The fluxes of the integrated continuum and line intensities 
(line plus continuum) within the big and small boxes are reported 
in Tab.[^ This corresponds to an observed 2 — 12 keV luminosity 
of L 2-12 = 3.4 X 10^® erg s“^ and L 2-12 = 2.6 x 10^® erg s“^ 
for the big and small boxes, respectively, assuming a distance of 
8 kpc to the Galactic center (Reid 1993; Reid et al. 2009). 


8.2 A new X-ray filament, XMM 0.173-0.413 

We observe a new X-ray filament extending ~ 2.2 arcmin per¬ 
pendicular to the Galactic plane and situated at I = 0.173°, 
b = —0.413°, which is almost directly toward negative longitudes 
from the GC Radio Arc. It coincides with the brightest segment 
of a much longer radio filament that extends toward the southern¬ 
most extensions of the filaments of the Radio Arc (Yusef-Zadeh et 


® To perform this task we use webpimms: httpsiheasarc.gsfc.nasa.govcgi- 
binToolsw3pimmsw3pimms.pl. As explained in § 2.1, the combined count 
rate is the sum of the EPIC-pn plus the EPIC-MOS count rates after scaling 
the latter exposure maps by 0.4. We then use the EPIC-pn, medium-filter, 
rate-to-fiux conversion computed with webpimms. 


al. 1989; 2004; see figures 7a, 16b,c, and 17a,b,c of the latter ref¬ 
erence, where the filament is labelled ”S5”), but because it is not 
continuous or exactly parallel with the filaments of the Arc, it is 
not completely evident that it is an extension of the Arc in three di¬ 
mensions. The X-ray filament has a hard X-ray colour and does not 
appear in the soft line images, indicating a non-thermal emission 
spectrum. 

Three other nonthermal radio filaments have been found 
to have X-ray emission along some portion of their lengths: 
G359.54-h0.18, G359.89-0.08, and G359.90-0.06 (Lu, Wang & 
Lang 2003; Sakano et al. 2003; Lu et al. 2008; Johnson et al. 2009; 
Morris, Zhao & Goss 2014; Zhang et al. 2014). XMM 0.173-0.413 
is the only one of the four known cases where the X-ray emission is 
not at or near a location where the radio filament shows unusually 
strong curvature. 

8.3 SNR excavated bubbles within the CMZ? 

The top and bottom panels of Fig.show the X-ray (2.5-4.5 keV) 
and 850 /im (Pierce-Price et al. 2000) maps of the Sgr B1 region. 
Four enhancements of X-ray emission are clearly present (shown 
by the green dashed ellipses in Fig. [^. In particular, G0.570- 
0.001, GO.52-0.046 and GO.40-0.02 correspond to holes in the dust 
distribution derived from the 850 /xm radiation (see bottom panel 
of Fig. [^. To reinforce this evidence, we observe that their X-ray 
edges can be traced in the dust distribution all around the X-ray en¬ 
hancements, suggesting a tight connection between the two. Such 
phenomenology is typical of SNe exploding within or near molecu¬ 
lar clouds and interacting with them, creating bubbles in the matter 
distribution (Ferreira & de Jager 2008; Lakicevic et al. 2014). In¬ 
deed, in this case, the SN ejecta might have cleared the entire region 
that is not filled with hot. X-ray emitting plasma, pushing away the 
ambient molecular material. However, this is not the only possibil¬ 
ity. The apparent X-ray enhancements could have resulted instead 
from a higher obscuration surrounding the submm holes, leading to 
higher X-ray extinction at the edges. We note that none of these re¬ 
gions has a known radio SNR counterpart (see Fig.[^. However, 
the radio emission might be confused within the very high radio 
background of diffuse emission created by G0.30-i-0.04, the several 
HII regions present in this region (see Fig. [^, and the bright, ex¬ 
tended synchrotron background of the GC. Enhanced X-ray emis¬ 
sion has already been reported towards GO.40-0.02 (Nobukawa et 
al. 2008) and close to GO.570-0.001 (see cyan ellipses in Fig.p^. 

8.3.1 Spectral analysis 

To further investigate the origin of these structures, we extracted a 
spectrum from each of these features (in either obsid 0694641301 
or 0694641201). We fitted each spectrum with a model composed 
of a SNR emission component (fitted with a PSHOCK model; a con¬ 
stant temperature, plane-parallel plasma shock model, meant to re¬ 
produce the X-ray emission from a supernova remnant in the Sedov 
phase) plus the emission components typical of the GC environ¬ 
ment such as a hot thermal plasma (with temperature in the range: 
kT = 6.5 — 10 keV), and an Fe Ka emission line, all absorbed 
by foreground neutral material (PHABS*(PSHOCK-i-APEC-i-GAUS) 
in XSpec). We assume that all these components have Solar abun¬ 
dances. Possible confusion effects produce uncertainties associated 
with the determination of the correct sizes of these candidate SNRs, 
so some of the results presented here, such as the dynamical time- 
scales, could thereby be affected. 
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2-4.5 keV emission 


GO.30+0.04 


00.224-0.032 




Galactic longitude 


Figure 17. (Top panel) 2-4.5 keV map. The dashed white ellipse shows the location of the radio SNR G0.30-i-0.04. The dashed and solid cyan ellipses show 
the positions of known X-ray SNRs and superbubbles, respectively. The green dashed ellipses indicate enhancements of soft X-ray emitting gas (GO.40-0.02 
was already observed in X-rays by Nobukawa et al. 2008), the magenta ellipse is used in the spectral analysis as a background region (Back in Tab. 7). The 
dashed yellow ellipse shows the location of the massive stars that could be in the tidal tail of the Quintuplet cluster (Habibi et al. 2013; 2014). (Bottom panel) 
850 pm map of the GC obtained with the SCUBA bolometer (Pierce-Price et al. 2000). The four dashed ellipses indicated by the green ellipses and showing 
enhanced X-ray emission correspond to holes in the 2mm gas distribution. 
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Big box 

Small box 

Flux 

Surf Bright. 

Flux 

Surf. Bright. 

7^1-2 keV = 19.0 

/ 1-2 keV = 6.2 

7^1-2 keV = 12.9 

fl-2 keV = 6.9 

F 2 - 4.5 keV = 155.0 

/2-4.5 keV = 50.7 

F 2 - 4.5 keV = 119.3 

/2-4.5 keV = 64.3 

F 4 . 5-12 keV = 290.7 

/ 4 . 5 -I 2 keV = 95.1 

F 4 . 5 -I 2 keV = 219.0 

/ 4 . 5 -I 2 keV = 118.0 

Esi xiii 4.4 

fSi xiii ~ 1-5 

Esi xiii — 3.4 

fsi xiii ~ 1-3 

F’Sa.. = 15.1 

fs XV = 4.9 

Es XV = 12.0 

fs XV = 6.4 

Eav xvii ~ 15.1 

f Ar xvii 4.9 

Eav xvii ~ 11.9 

f Ar xvii ~ 6.4 

Eca xix — 18.4 

fCa xix — 9.9 

Eca xix — 14.3 

fCa xix 7.7 


Ta ble 6. Fluxes (F) and surface brightnesses (f) of the continuum and the line intensities (line plus continuum) integrated over the small and big boxes described 
in § |8.l| The fluxes are given in units of 10“^^ erg cm“^ s“^, while the surface brightnesses are given in 10“^^ erg cm“^ s“^ arcrnin”^. 


GO.40-0.02 shows a best-fit temperature and normalisation of 
the warm plasma associated with the SNR component of kT — 
0.55 zb 0.1 keV and ^psho = 2.4^^ x 10“^, while the column 
density is observed to be Nh — 7.7 zb 0.8 x 10^^ cm“^ (see 
also Nobukawa et al. 2008). To test whether the observed X-ray 
enhancement is due to a real variation of the intensity of the soft 
X-ray emission or whether it is the product of lower extinction, 
we also extracted a spectrum from a nearby background compar¬ 
ison region (magenta in Fig. and ’’Back” in Tab. |^. This sec¬ 
ond region has the same size as GO.40-0.02 and it is located in 
a fainter region in X-rays, characterised by higher Nh, as sug¬ 
gested by the 850 /xm map (see bottom panel of Fig. [^. This 
background region shows a slightly higher absorption column den¬ 
sity, Nh ~ 8.8 zb 1 X 10^^ cm“^, and no significant warm plasma 
component. If we impose the presence of a warm plasma compo¬ 
nent having the same temperature and r (the ionisation timescale 
of the shock plasma model) as observed in GO.40-0.02 we obtain 
an upper limit to its normalisation of Aphsho < 6x10“^. This sug¬ 
gests that the enhanced X-ray emission towards GO.40-0.02 is due 
to a real excess of X-ray emission and is not a simple byproduct of 
lower extinction. The thermal energy, the dynamical timescale and 
the size of GO.40-0.02 are Eth ~ 1.9 x 10^° erg, tdy ~ 3700 yr 
and 8.6 X 5.5 pc^, respectively, as expected for a young SNR in the 
Sedov-Taylor phase (derived from the equations shown in Maggi et 
al. 2012). 

Similar parameters characterise GO.52-0.046 (kT — 0.77 zb 
0.3 keV, Apsho = X 10■^ Nh = 7.9 ± 1.1 x lO^^ cm-^). 
Therefore, this feature also appears to be consistent with an SNR 
origin. However, although we derive a dynamical age of the same 
order (tdy ~ 1700 yr), the energy inferred for the SN explosion is 
substantially lower, Eth ~ 5 x 10^^ erg. 

The spectrum of GO.57-0.001 is characterised by significantly 
lower statistics and a higher column density of absorbing mate¬ 
rial. The best fit prefers a low temperature plasma, at the limit 
of detection. We fix its temperature to a relatively low value of 
kT = 0.6 keV and find Apsho = 4.31“ x 10“^ and Nh = 
9.5 zb 2 X 10^^ cm“^. The derived thermal energy and dynamical 
times are Eth ~ 2.6 x 10^^ erg and tdy ~ 1600 yr, respectively. 
The soft X-ray excess of this feature lies spatially very close to a 
region of X-ray excess that has been traced by Fe XXV line emis¬ 
sion (Nobukawa et al. 2008). We also note that the region defining 
GO.57-0.001 almost completely contains a diffuse X-ray source de¬ 
tected both by the Chandra and A^CA satellites (Senda et al. 2002). 
The Chandra image shows a very compact (~ 10” radius) and hot 
shell, GO.57-0.018, possibly the youngest SNR in the Galaxy, less 
than about 100 yr old. However, Renaud et al. (2006), because they 


found neither radio nor nucleosynthetic decay products (such as 
^^Ti), questioned such an interpretation. Further investigations are 
required to understand the link, if any, between these features. 

The morphology of GO.224-0.032 appears more complex, 
compared to the other SNR candidates. The edge of the X-ray emis¬ 
sion is well defined only towards the Brick molecular cloud (des¬ 
ignated MO.25-1-0.01 in Fig. 5) that, with its very high column den¬ 
sity, can obscure the soft X-ray emission there. In any case, the 
fit of its X-ray spectrum shows parameters typical of a SNR. In 
particular, we derive a thermal energy and a dynamical time of 
Eth ~ 2.6 X 10^° erg and tdy ~ 1800 yr, respectively. There¬ 
fore, GO.224-0.032 might be a new SNR partly obscured by the 
brick molecular cloud. In such a case, the true size and the energy 
estimate are likely larger. 

Overall, we remark that, if these SNR candidates are real, their 
dynamical timescales are extremely short, which would imply an 
extremely high supernova rate. The SN rate in the CMZ has been 
estimated to he as high as 0.4 SN per millenium (Crocker et al. 
2011). However, we caution that our dynamical time-scales could 
be off either because of a higher ambient density than we have as¬ 
sumed, or because absorption or confusion effects do not allow us 
to distinguish the proper border of the SNRs, or to detect any colder 
and more extended portions that might be present. 


8.3.2 Expanding molecular shells 

The observed temperatures, ages and sizes of these SNR candidates 
are consistent with a Sedov-Taylor framework expanding into an 
average ambient density between 1 — 10 cm“^ (higher density en¬ 
vironments would result in older and cooler SNRs; Ostriker & Mc¬ 
Kee 1988). 

If the X-ray enhancements described above truly arise from 
SNRs interacting with and carving bubbles inside or near the sur¬ 
faces of molecular clouds (preceded, perhaps, by the wind of the 
massive progenitor), we should observe clear traces of such events 
also in the kinematics of the surrounding molecular matter. Such 
a complex and delicate investigation is beyond the scope of the 
present paper. Nevertheless, we note that Tanaka et al. (2009) dis¬ 
covered an expanding SiO shell (SiOO.56-0.01) centered at I ~ 
0.56°, b ~ —0.01° and having a size of ~ 3.0 x 3.4 pc^. The 
center and size of the expanding SiO shell closely match the peak 
and size of the X-ray emission of GO.57-0.001 and suggest an as¬ 
sociation between the two. In particular, high-velocity clumps have 
been found consistent with the idea that the SiO shell consists of 
swept-up material. Tanaka et al. (2009) calculated a kinetic energy 
of Ekin ~ 10^°“ erg for SiOO.56-0.01. This strongly suggests that 
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parameter 

unit 

GO.40-0.02 

GO.52-0.046 

GO.570-0.001 

GO.224-0.032 

Back 

kT 

keV 

0.55 ±0.1 

o 

1 + 

3 O 

to ^ 

0.6t 

0.54 ±0.1 

0.55t 

■^psho 


2.4+^ X 10-2 

5^3 X10-3 

4.3iy X 10-3 

3^2 X10-2 

< 6 X 10-3 

rtt 

s cm“^ 

> 1.7 X lO^i 

> 7.5 X 10^0 

> 1.4 X IQio 

> 3.5 X lO^i 

7.5 X 10“t 

Size 

pc 

8.6 X 5.5 

5.9 X 2.7 

4.2 X 2.1 

5.4 X 2.6 

8.6 X 5.5 

Nh 

10^^ cm“^ 

7.7 ±0.8 

7.9 ±1.1 

9.5 ±2 

7.4 ± 1 

8.8 ± 1 

X^/dof 


807/777 

342/364 

141/134 

369/338 

834/741 

rie 

cm“^ 

1.4 

1.4 

2.1 

11 


i'dy 

yr 

3.7 X 10^ 

1.7 X 103 

1.6 X 103 

1.8 X 103 


Eth 

erg 

1.9 X 10^0 

5.0 X 10^9 

2.6 X 10^9 

2.6 X 10^0 



Table 7. Best-fit and derived parameters of the SNR candidates described in § |8.3| fParameter unconstrained. JValue weakly constrained by the high column 
density of neutral material, therefore fixed for the corresponding fit. ft Ionisation time-scale of the shock plasma model. 


GO.57-0.001 is indeed a SNR caught in the process of carving its 
bubble. Further studies of the gas kinematics around the other X- 
ray enhancements are required to establish their real nature. 

We note that the thermal energy estimated for these SNRs is 
observed to be systematically lower than the theoretical value for 
the remnant of a standard type II SN expanding into the interstellar 
medium. This might result from a relatively higher ambient density 
in the GC, leading to greater energy dissipation, or from a signif¬ 
icant fraction of the energy budget going into the inflation of the 
bubbles and the production of cosmic rays. 

We also note that GO.570-0.001, GO.52-0.046 and GO.40-0.02 
are located within the trail of massive stars that have been hypoth¬ 
esised to have tidally escaped from the Quintuplet cluster (see Fig. 
[^and Habibi et al. 2013; 2014). This raises the possibility that 
some of these SNRs might be associated with SN explosions from 
stars originating in this massive, young stellar cluster. 

8.4 Origin of the Sgr A X-ray lobes 

All the soft X-ray maps (see Figs. [TT] and show 
the presence of two extended features, with a size of roughly 5 — 
10 pc, located to the Galactic north and south of Sgr A*, the so- 
called ’’bipolar Sgr A lobes” (Morris et al. 2003; 2004; Markoff et 
al. 2010; Heard & Warwick 2013). 

8.4.1 Lobe morphology 

The lobes appear to have roughly oval shapes with co-aligned ma¬ 
jor axes oriented perpendicular to the Galactic plane. They ap¬ 
pear joined at the position of Sgr A* suggesting the latter is their 
point of origin (see Fig. |^. The top panel of Fig. [II]andFig.|7] 
show that the lobes’ emission is characterised by a smaller ratio of 
soft X-ray lines to continuum (therefore characterised by a greener 
colour) compared to the surrounding regions (appearing with a 
redder colour) such as the superbubble, GO. 1-0.1 and the Radio 
Arc (Fig.[TT][^and[^. This suggests that the lobes, although they 
show thermal emission lines (see §4.2| ), have either a stronger non- 
thermal component or smniflcantly hotter thermal emission than 
the surrounding regionq^ We also note that the eastern portion of 
what appears to be part of the southern lobe has a colour as red as 

A hot plasma, with temperatures of ~ 2 — 4 keV, produces intense X- 
ray emission but weaker soft X-ray lines, compared to a plasma having a 
temperature around 1 keV. 


GO. 1-0.1 and the superbubble regions (see Fig. 10). Therefore, this 
emission might not be associated with the lobes, but rather with 
GO. 1-0.1 or the edge of the superbubble (however, a gap such as 
might be produced by a foreground dust lane, appears to separate 
the lobes’ emission from GO. 1-0.1). Figures |7] andalso show that 
the surface brightness of the northern lobe decreases with distance 
from Sgr A* (Heard & Warwick 2013). 

The bottom panel of Fig. pT] shows that the lobes have an or¬ 
ange colour, indicating harder soft X-ray emission compared to the 
surrounding regions. In particular, a brighter and harder (yellow- 
green) linear structure outlines the northern lobe, shaping it to have 
a well-deflned and symmetric spade structure, with a sharp transi¬ 
tion at the border. The sharpness of the transition suggests the pres¬ 
ence of a limb-brightened shock, indicating that the lobe is a bubble 
enclosed by a thin shell of hot, compressed material. This claim is 
strengthened by our analysis of the Chandra data (see Fig.[^and[^ 
Baganoff et al. 2003; Lu et al. 2008; Muno et al. 2008). The supe¬ 
rior Chandra spatial resolution, in fact, allows us to note that these 
projected linear features are running right along the lobes’ edges, 
indeed confirming the presence of a shock (Fig. |^and|^. We also 
note that the emission from the northern half of the lobe seems to 
be mainly due to three harder filaments converging at the top in a 
cusp, having radio continuum counterparts (Zhao et al. 2015) and 
associated Paschen-a emission, indicating that these are thermal 
features (see Fig.p^andp^. 

In the southern lobe, two bright knots are observed in the cen¬ 
ter and at the tip. Interestingly these appear to be located approxi¬ 
mately at the same distance and in the opposite direction, compared 
to Sgr A*, as two enhancements present in the northern lobe. More¬ 
over, the two bright knots have a green-yellow colour (upper panel 
of Fig. 11) similar to their apparent counterparts in the northern 
lobe. This suggests both: i) a similar physical origin for these fea¬ 
tures in both the north and south lobes and; ii) that the process that 
created the lobes is symmetric about the Galactic plane and its en¬ 
gine is (or is located close to) Sgr A*. The obvious interpretation 
of this morphology is that energetic events simultaneously ejected 
diametrically opposed blobs of hot gas. 

However, upon closer inspection, the northern and southern 
lobes do not appear completely symmetric. For example, compared 
to the northern lobe, the western side of the southern lobe and the 
region close to Sgr A* appear suppressed (see Fig.[^|^and discus¬ 
sion following). On the other hand, the eastern side of the south¬ 
ern lobe appears to extend further (e.g. further east compared to 
G359.977-0.076) than the corresponding boundary of the northern 
lobe (located close to e.g. G359.974-0.000). As described before. 
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the emission around and east of G359.977-0.076 has a different 
colour and might therefore be associated with either the superbub¬ 
ble, with GO. 1-0.1, or it could be a feature that is independent of 
either of these and of the lobes. 

The bottom panel of Fig. shows that the region with de¬ 
pressed soft X-ray emission south of Sgr and on the western 
side of the southern lobe spatially coincides with the presence of 
the 20 km s“^ molecular cloud, which is thought to be located in 
front of Sgr A* (Coil et al. 2000; Ferriere 2009) and to have a large 
column density. Soft X-ray emission could be produced there but 
be completely obscured to us by this intervening cloud (see Fig.[^. 
To reinforce this idea, we note that, in fact, at this position, hard 
X-ray radiation (between 4.5 and 12 keV) is observed by Chan¬ 
dra to have a non-thermal spectrum and to be extended (Morris 
et al. 2003). Therefore, this hard radiation, which is able to pene¬ 
trate the cloud, could be produced by strong shocks at the bubble’s 
border. Similar hard non-thermal filaments are observed in several 
places at the border of the northern lobe (Morris et al. 2003). This 
strongly suggests that the actual border of the southern lobe is lo¬ 
cated further west than the images reveal, and that the lobe’s soft 
X-ray emission is obscured there. Therefore, once the effect of ab¬ 
sorption by molecular clouds (e.g. the 20 km s“^ cloud) is con¬ 
sidered, the symmetry between the northern and southern lobes 
appears more clearly. The linear or filamentary structures (such 
as G359.974-0.000, G359.970-0.009, G359.959-0.027, G359.945- 
0.044, G359.942-0.045, 359.933-0.039, see Fig.|8} observed in the 
northern lobe might be present in the southern one as well, but be 
suppressed by the intervening absorption. 

Figure shows the comparison between the Chandra and 
Paa emission. The Paa map clearly show the presence of tendrils 
of foreground absorbing material, running north-south along the 
extension of the lobe. Interestingly, the same tendrils are also evi¬ 
dent as absorption lanes in the X-ray image. We also note that two 
Paa emitting luminous stars are contributing to, if not dominat¬ 
ing, the ionisation of the gas surrounding them. However, there is 
a close correspondence between the Paa emission and the soft X- 
ray emission at the northern tip of the lobe, so strong shocks might 
be contributing both to the ionization and to the production of very 
hot post-shock gas that can emit X-rays. Indeed, we hypothesize 
that some portion of the hot wind that we see in the softest X-ray 
bands, presumably emanating from near Sgr A*, is undergoing a 
shock where it encounters ambient interstellar material at b ~ 0.02 
degrees, and that it is thereby blocked from continuing to higher 
latitudes. That shock is manifested both as a horizontal feature in 
the Paa image, and as a diminution in the brightness of the X-ray 
emission proceeding north from that latitude. There is still some 
weaker, extended X-ray emission at higher latitudes, indicating that 
not all of the outfiowing wind is completely blocked. In addition, 
the portion of the X-ray emitting plasma lying behind the shock 
front appears to be absorbed along the shock front, creating a hori¬ 
zontal shadow in the extended X-ray emission, and suggesting that 
the shock at b ~ 0.02 degrees has created a thick, compressed layer 
that absorbs the X-rays coming from behind it. 

8.4.2 Lobes collimated by the Circumnuclear Disk (CND) 

Morris et al. (2003) observed that the CND has a size and orien¬ 
tation that are consistent with it being the agent that collimates an 
isotropic outflow from the Sgr A* region, thereby creating the bipo¬ 
lar lobes of hot plasma. Those authors further suggested that the 
sequence of enhancements along the axis of the lobes might have 
resulted from a series of energetic mass ejections from the imme¬ 


diate environment of Sgr A*. A similar scenario has recently been 
discussed by Heard & Warwick (2013b). At an outflow velocity of 
10^ km s“^, it would have taken 2 — 6 x 10^ yr to infiate the lobes. 
Assuming a thermal emission model, an electron density Ue in the 
range 1 — 10 cm“^ can be inferred (Morris et al. 2003; Heard & 
Warwick 2013b) and, assuming that the line-of-sight depth is equal 
to the projected width, the total hot plasma mass involved in the 
X-ray lobes is only about 1 — 3 Solar masses. 

8.4.3 Energetics of the Lobes 

Except for the broad intensity enhancements along the axis of the 
lobes, the surface brightness is relatively smoothly distributed. As¬ 
suming a continuous and constant outflow, we examine the energy 
budget of the lobes. In particular, integrating the measured energy 
density over a cylinder of 5 pc radius and 12 pc height (the approx¬ 
imate sizes of the lobes), a thermal energy of Eth ~ 9 X 10^® erg 
is estimated. 

It has been estimated that the massive stars in the central par¬ 
sec collectively lose ~ 5 x 10“^ M© yr“^ in stellar winds, with 
velocities ranging from ~ 300 to 1000 km s“^ (Geballe et al. 1987; 
Najarro et al. 1997; Paumard et al. 2001). The total kinetic energy 
thermalised by shocks is then ~ 5 x 10^® erg s“^ for such a mass 
outfiow rate and an outfiow velocity of 1000 km s“^ (Quataert & 
Loeb 2005). Therefore the energy released within the time needed 
to infiate the lobes (~ 4 x 10^ yr) is equivalent to ~ 5 x 10^^ erg, 
therefore giving an important contribution to the generation of the 
lobes. 

The lobes might also be traceable to the accretion fiow onto 
Sgr A*. We note that, as estimated by Wang et al. (2013), only 
~ 1 % of the matter initially captured at the Bondi radius presently 
reaches the innermost regions around Sgr A*. The rest of the accre¬ 
tion power, estimated to be ~ 10^^ erg s“^ (Wang et al. 2013), is 
probably converted to kinetic energy and used to drive an outfiow 
that carries away the bulk of the infiow, sculpting the environment 
with its ram pressure. If so, within the lobe inflation time, a total 
energy of ~ 10^° erg would have been deposited. 

Within the CMZ, X-ray refiection nebulae indicate that a few 
hundred years ago Sgr A* was more active, being ~ 10® times 
brighter than at present for approximately 5 — 10 % of the time in 
the past millennium (experiencing Lx ~ 10®® erg s“^; see Ponti 
et al. 2013 for a review). Could the lobes have been created by sim¬ 
ilar events that occurred over the past 10^ yr? The light crossing 
time of the CMZ limits our capability to trace Sgr A*’s past activ¬ 
ity beyond about 10® yr ago, so it is difficult to directly trace the 
echoes of possible energetic events on such time scales. However, 
if the process has been active over the past (5 — 10)xl0®yrat 
roughly the same rate, (therefore active at F ~ 2 x 10®® erg s-^ 
for ~ 10® yr in the past 10^ yr), a total integrated energy equal to 
~ 5 X10^® erg should have been generated. If Sgr A* ’s past activity 
was characterised by outbursts with associated outflows having ki¬ 
netic luminosity comparable to the radiated power (therefore much 
higher than in soft state stellar mass black holes; Ponti et al. 2012), 
then these events could be the primary source (or at least contribute) 
to form the lobes. 

All these processes appear similarly likely to have an impact 
on the formation of the lobes, from an energetic point of view. 
However, we note that the first two mechanisms are powered by 
a quasi-continuous outfiow from Sgr A* or from the central stel¬ 
lar cluster. In such scenarios, therefore, the sharpness of the edges 
at the extremities of the lobes remains rather puzzling, favouring 
explosive-outbursting scenarios. 
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Figure 18. (Left panel) Chandra RGB image of the northern lobe. See Figure 8 for more details. (Right panel) Pa-a image of the northern lobe, tracing 
thermal, ionized gas (from Wang et al. 2010). Continuum emission from stars has been removed (Dong et al. 2011), so the only stars that appear are those that 
have strong Pa-a emission lines. Two Paa-emitting luminous stars located at I ~ 359.935°, b ~ 0.21° and I ~ 359.925°, b ~ 0.45° (Mauerhan et al. 2010; 
Dong et al. 2012) are probably responsible for at least part of the ionisation indicated by the Pa-a. See www.mpe.mpg.de/heg/gc/ for a higher resolution 
version of these figures. 


8.4.4 Are SgrA^s lobes the SNR ofSGR J1745-2900? 

The recent discovery of a young magnetar, SGR J1745-2900 (De- 
genaar et al. 2013; Dwelly & Ponti 2013; Mori et al. 2013), most 
probably in orbit around Sgr A* (Rea et al. 2013), raises the quest 
for finding its young SNR. SGR J1745-2900 is estimated to be only 
about 9 X 10^ yr old and to be located at <0.07 — 2 pc from 
Sgr A* (Rea et al. 2013). Therefore, the supernova that generated 
SGR J1745-2900 should have exploded near the centroid of the 
lobes and likely inside the inner radius of the CND. If the shock 
propagated at the sound speed (Vs ~ 750 km s“^), then a present 
size of ~ 6.5 pc would be expected. This is slightly smaller than 
the observed size of the lobes (~ 12 pc), but is consistent with 
the observed size within the uncertainties in the age and the sound 
speed. We also note that the thermal energy content in the lobes 
(Eth ~ 9 X 10^^ erg) is lower, but of comparable order of mag¬ 
nitude, to the energy released by typical supernova explosions. It 
is therefore plausible that the lobes are indeed the SNR associated 
with SGR J1745-2900. 

Another viable possibility is that the lobes have been gener¬ 
ated by the supernova that created the PWN candidate G359.945- 
0.044, located only ~ 8” from Sgr A*, with an estimated age of 
few thousand years (Wang et al. 2006). It is therefore possible that 


at least one SN exploded close to Sgr A* within the last ~ 10 kyr. If 
so, its blast wave likely propagated into a pre-existing hot, low den¬ 
sity cavity created by Sgr A*’s outfiows and the collective winds of 
the central stellar cluster. Given the density and the temperature in 
the lobes (kT ~ 2 keV; Morris et al. 2003), the shock is estimated 
to reach 15 pc in 9 x 10^ yr (Wang et al. 2005), a value consistent 
with the observed size of the lobes. 

As mentioned above, the presence of the sharp edges to the 
lobes seems to favour an explosive mechanism for their creation. 
A supernova exploding inside the CND would expand into the pre¬ 
existing, stationary outflow from the center and be collimated in 
the same way. Hydrodynamical simulations typically show that SN 
shock fronts are reflected away when encountering the walls of 
a dense molecular cloud, such as the CND (Ferreira & de Jager 
2008). In this scenario the sharp edge of the lobes would be due to 
the SN shock front. If the supernova recurrence time is longer than 
or comparable to the lobe expansion time (a few thousand years) 
then this would not appear as a stationary process. Assuming a re¬ 
currence time between 1 —10 x 10^ yr (similar to the SN recurrence 
time of the central young cluster, in which ~ 100 massive stars pre¬ 
sumably becoming SN over a ~ 10^ yr time interval), we estimate 
a time-averaged kinetic power release of 3 — 30 x 10“ erg s“^ This 
indeed suggests that: i) SN explosions of the central star cluster can 


© 0000 RAS, MNRAS 000, [T]-?? 



29 


contribute to powering the lobes; ii) the lobes are quasi-stationary 
features; and iii) it is not unlikely that we observe such features 
created by a rare event such as a SN explosion. 

Finally, we note that, although the characteristics of the X- 
ray emission from the lobes appear consistent with being the X-ray 
remnant of a recent supernova that exploded within a few parsecs 
of Sgr the lack of associated nonthermal radio emission from 
such a young SNR is problematical for this hypothesis. 

8.5 G359.77-0.09 and G359.79-026: a ring from a hot 
superbubble southwest of Sgr A* 

A series of diffuse, soft features appear to the southwest of Sgr A* 
(see Fig. [^, namely G359.77-0.09, G359.79-0.26 and a newly 
recognised extended feature, GO.0-0.16. If these are physically con¬ 
nected, they form, in projection, a roughly elliptical shape whose 
major axis has an inclination of about 60° with respect to the Galac¬ 
tic plane (see Fig.[^[6l[7| p^pT]p^ andpl). These features show 
similar colours and strong, soft line emissions, indicating a simi¬ 
lar thermal origin (see Fig.[^[^and[TTJ. This elliptical structure 
appears brightest at the softest energies, however it is not observed 
below 1.5 keV, suggesting a location near the Galactic centre and a 
low temperature for its plasma, compared to the surrounding emis¬ 
sion. We note that this feature is characterised by a very bright edge 
with strong Si XIII emission on the outside of the ellipse, suggest¬ 
ing a lower temperature of the edge compared to the interior. 

The ellipse center is located at I ~ 359.9°, b ~ —0.125° and 
it has minor and major axes of about 7.8 and ~ 12 arcmin, respec¬ 
tively (corresponding to 18 and 28 pc). Both Mori et al. (2009) and 
Heard & Warwick (2013) considered that these structures/group 
of structures were physically connected and form a superbubble 
candidate. We note that the recognition of such an elliptical ring 
critically depends on the presence of the dark lane running from 

I ~ 0.02°, b -0.22° to I ~ 0.05°, b -0.07° (see Fig.|^ 

[^[TT]|7] and [B, which separates GO.0-0.16 from the emission 
of GO. 1-0.1; this dark lane helps define the quasi-continuous ellip¬ 
tical morphology of the ring. However, the lane might simply be 
due to absorption by foreground material, in which case GO.0-0.16, 
forming the eastern part of the ring, could simply be connected to 
GO. 1-0.1. 


8.5.1 S XV emission filling the superbubble 

The S XV emission provides a key piece of information to better un¬ 
derstand the superbubble. We note, in fact, that the S XV emission 
completely fills the region inside the ring with a roughly uniform 
brightness (several times brighter than in the surrounding region) 
and sharply drops at the ring’s edge. This indicates that the super¬ 
bubble is a shell of hot gas that we see projected onto the plane 
of the sky (see middle top panel of Fig. [^. To further corrobo¬ 
rate this, we note that a sharp emissivity drop appears to be located 
just outside of the ring, running all around the ring’s external edge. 
Such X-ray depression might be produced by a high column den¬ 
sity of cold gas pushed away by the superbubble’s shock front and 
accumulated in large quantities just outside the shock. If so, the 
observed depression could be indicating that the superbubble is lo¬ 
cated in front of GO. 1-0.1. Fig. also shows a small depression in 
the top part of the northern lobe that could easily be explained by 
absorption associated with the superbubble if it is located in front 
of the lobes. This situation would then be somewhat analogous to 
the colour variation in the south lobe. 


Mori et al. (2009) and Heard & Warwick (2013) estimate a 
total thermal energy contained within the superbubble of Eth ~ 
10^1 ej-g (where / is the volume filling factor of the emitting 
plasma). Such a large energy content does, indeed, require multiple 
supernova events. Those authors also estimate for G359.77-0.09 
and G359.79-0.26 an ionisation time-scale of tion ~ 3 x 10^ yr 
(assuming / ~ 1). 

8.5.2 Origin 

The origin of such a superbubble is not clear. We note that many 
of the massive stars that are suggested to have escaped from the 
Quintuplet cluster (Habibi et al. 2013; 2014) are projected inside 
the superbubble. It is therefore possible that explosions of stars lost 
by the Quintuplet cluster have contributed to energising the super¬ 
bubble. 

A more speculative point is that the estimated age of the su¬ 
perbubble is of the same order of magnitude as the recurrence time 
of tidal disruption events by Sgr A^: Itde ~ 1 — 3 x 10^ yr 
(Alexander & Hopman 2003). While Sgr A* appears, in projec¬ 
tion, to be located inside the superbubble, it is ~ 6.7 arcmin off 
from the superbubble’s center. Khokhlov & Melia (1996) suggested 
that an explosion associated with a tidal disruption event would 
liberate a large amount of energy on the order of E ~ 10^^ erg 
that would propagate as a powerful shock wave into the local in¬ 
terstellar medium. As with the remnant of the SGR J1745-2900, 
we expect that the shockwaves of a tidal disruption event would 
interact with the CND. However, the unbound part of the tidally 
disrupted star would be ejected into a limited solid angle, produc¬ 
ing a strongly elongated and asymmetric remnant (Khokhlov et al. 
1996; Ayal et al. 2000). Such a remnant would then appear as a 
very energetic shell of hot gas and remain visible for a time compa¬ 
rable to the age of a typical SNR. Assuming a shock survival time 
of- 1-10 X 10^ yr, we could potentially observe a few remnants 
resulting from tidal disruptions. We suggest that the superbubble 
G359.77-0.09 has properties that make it a possible candidate. No 
other feature with properties obviously related to a tidal disruption 
event appear to be observed close to Sgr A*. 

8.6 The arc bubble: a second superbubble in the GC 

Highly enhanced soft X-ray emission is observed east of Sgr A* 
from the region called GO. 1-0.1. This feature appears in Figures 
and m as a slightly elliptical feature of enhanced emission with 
center at / — 0.09°, 5 — —0.09° and with radius of — 5 arcmin 
(corresponding to — 10 pc). The top panel of Fig. m shows that 
GO. 1-0.1 and the Radio Arc regions both show distinct red emis¬ 
sion. This indicates large equivalent widths of the emission lines 
from this plasma and therefore a thermal origin. However, the top 
panel of Fig.[^and the bottom of Fig.pTjshow strong colour gradi¬ 
ents within these regions, indicating that they might have different 
contributions from distinct components. 

The PWN candidate GO.13-0.11 (/ = 0.131°, b = -0.111°; 
Wang et al. 2002; Heard & Warwick 2013) stands out from the gen¬ 
eral thermal emission in GO. 1-0.1, appearing as a distinct point-like 
source characterised by a light blue/white colouf^ indicating its 
non-thermal origin (dominated by intense soft and hard continuum 

Note that in the top panel of Fig. [^despite an enhancement of diffuse 
emission from the region surrounding the core of GO. 13-0.11, no point like 
source is detected in the soft X-ray line image. 
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emission with no soft line emission, see the top panel of Fig. ll|^ 


The point-like head of GO. 13-0.11 appears to be accompanied by a 
tail extending to the south for 4.5 — 5 pc; in Fig.[^ it appears with 
a white-violet colour. 

Heard & Warwick (2013) suggest that the SN that generated 
GO. 13-0.11 might be the source of the soft X-ray emission from this 
region. Those authors present a spectral study of the X-ray emission 
from GO. 1-0.1 and find a gas temperature of kT = 1.1 zb 0.1 keV, 
and a column density of Nh = 5.6 zb 0.5 x 10^^ cm“^, indicat¬ 
ing a GC location of this emission, and abundances that are about 
1.8 times Solar. Assuming that the plasma volume is only 3.5 pc^, 
corresponding to only 1.5 arcmin radius around the PWN (see the 
red circle in Fig. [^, the authors estimated a thermal energy of 
Eth = 3.1 X lO^^rg (and a plasma ionisation time-scale of at 
least t = 1.8 X 10^ yr), thus consistent with being produced by a 
single supernova explosion. 


8.6.1 Energetics of the arc bubble 

We note from Figs.p^[3l p^pT] [7]p^andpl]that GO. 1-0.1 extends 
further from GO. 13-0.11 than the 1.5 arcmin region size considered 
by Heard & Warwick (2013), with no clear boundary at 1.5 arcmin 
(see Fig.|^. To illustrate this, Fig.p^shows the soft emission lines 
RGB image with colour scales chosen to highlight intensity varia¬ 
tions present within this region. The left panel of Fig. shows the 
contours of the S XV emission overlaid on the 20 fim MSX image 
(Price et al. 2001). Figure [^clearly shows that the empty mid-IR 
bubble (the so called arc bubble; Levine et al. 1999; Rodriguez- 
Fernandez et al. 2001; Simpson et al. 2007) is completely filled 
with warm X-ray emitting plasma and that the soft X-ray emis¬ 
sion is not confined to within ~ 1.5 arcmin of GO. 13-0.11, but it 
extends much further, for about 7 arcmin. Assuming a uniform sur¬ 
face brightness, if the bubble is 4.5 times larger, we would expect 
a thermal energy of Eth ~ 1-5 x 10^^ erg, thus most probably 
requiring multiple supernova events, and supporting the notion that 
GO. 1-0.1 is a second superbubble candidate in the GC. 

Rodriguez Fernandez et al. (2001) have noticed that the radio 
Arc bubble is filled with continuum X-ray emission seen by ASCA 
which they ascribed to X-ray sources inside the bubble. Here we 
find that the bubble is in fact filled with diffuse thermal X-rays, 
most likely originating from SN explosions of massive stars asso¬ 
ciated with the Quintuplet cluster. 

We also note that the soft line emission is at least as extended 
as the arc bubble and is highly inhomogeneous (see Figs. 19 and 
20). Three depressions having roughly circular shapes can be dis¬ 
cerned in GO. 1-0.1 (Fig.|^. Two cavities are located at about the 
same latitude, with centers close to I = 0.057°, b = —0.067° 
and to I = 0.116°, b = —0.071° and with radii of ~ 1.6 arcmin 
(corresponding to ~ 3.7 pc) and ~ 1 arcmin, respectively. These 
cavities appear to be surrounded by a thin rim of brighter emission. 
A third depression is centered at I = 0.083°, b = —0.123°, with 
~ 1.8 arcmin radius. This cavity also seems to be confined by a 
thin shell of brighter material, except for its southern edge, where 
it appears open (see Fig.[^, possibly because of the presence of a 
dark absorbing lane. 


Please note that an unrelated point source atZ = 0.142°, 6 = —0.109° 
is located at the same latitude, but ~ 1.5 pc to the Galactic east of the PWN 
candidate GO. 13-0.11. 


8.6.2 Association with the Quintuplet cluster? 

Despite its offset position from the centre of the mid-IR arc bub¬ 
ble, the Quintuplet cluster is often considered responsible for cre¬ 
ating and maintaining the bubble with a combination of supernovae 
and strong stellar winds (e.g., Simpson et al. 2007). Johnson et al. 
(2007) invoke a possible non-uniformity of the ambient medium as 
a possible origin of this asymmetry. This might also apply to the X- 
ray emission. However, we note that the Quintuplet cluster is mov¬ 
ing supersonically within the CMZ (Stolte et al. 2014). Given its 
projected velocity, the cluster would have taken ~ 100 kyr to cross 
the width of the IR arc bubble, in which case the bubble would have 
been inflated on a time scale much smaller than typical superbubble 
formation times (Castor et al. 1975; Weaver et al. 1977; Me Low 
& McCray 1988). Several SN explosions in that amount of time 
would therefore have been required. Furthermore, the Quintuplet 
cluster would have been located in the middle of the two northern 
cavities about 4 x 10^ and 9 x 10^ yr ago, respectively (the right 
panel of Fig. [20|illustrates the direction of the cluster’s motion and 
the position of the cavities). The relatively small cavities observed 
in the northern part of GO. 1-0.1 are unlikely to have been generated 
by multiple cluster stars. In fact, a supernova exploding in the hot 
plasma of GO. 1-0.1 is expected to undergo a significantly differ¬ 
ent evolution than a typical SNR. In particular, the sound velocity 
is significantly larger than in a typical low pressure medium. Tang 
& Wang (2005) have shown that the shock velocity follows a Se¬ 
dov solution but quickly deviates from it when it becomes mildly 
supersonic. This translates into a much faster evolution and much 
larger cavities would be expected if the SN exploded 4x10^ and 
9 X 10^ yr ago, when the quintuplet cluster was at that location. It is 
more likely, therefore, that the two cavities to the Galactic west of 
the Quintuplet might have been generated by supernova explosions 
of massive stars either stripped from the Quintuplet cluster (Habibi 
et al. 2013) or having no association with it. 

The large thermal energy filling the arc superbubble could 
have been produced by some combination of winds from the young 
stars and by multiple supernova explosions, including the super¬ 
nova explosion associated with the PWN G0.13-0.11p^ 


8.7 A hot atmosphere around the GC - A link to the GC 
lobe? And to the Fermi bubbles? 

8.7.1 General morphology 

As observed in all soft X-ray maps (Figs. and pT]) and con¬ 
firmed by the soft plasma intensity map (obtained through spectral- 
images decomposition, Fig.[^, the regions at higher Galactic lati¬ 
tudes are significantly brighter in soft X-rays than the regions closer 
to the disc, presumably in part because of the smaller extinction at 
the higher latitudes. 

The western border of this enhanced emission is defined by 
a relatively sharp edge between I — 359.63°, b = 0.06° and 
I = 359.55°, b = 0.46° (Fig.[^ see also Figs.and [^. The 
soft X-ray emission peaks above the location of the GC Radio Arc, 
appearing as a continuation of the Radio Arc itself. Further west the 
soft X-ray emission appears to fade with increasing Galactic longi¬ 
tude. In particular, the spectral decomposition provides hints for the 
presence of an eastern edge, fainter but similar to the western edge, 
of this hot GC atmosphere. However the presence of an edge to 


We note that the PWN GO. 13-0.11 is located right in the middle of the 
mid-IR arc bubble. 
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Figure 19. RGB image with colour scale chosen to highlight enhancements and depressions in the diffuse emission east of Sgr A*. In red the sum of the Si-S 
and S XV bands is shown. In green the sum of S-Ar plus the S xv and Ar xvii bands is shown. The blue image shows the sum of the Ar-Ca plus Blue-Ca and 
Ca XIX bands (see Tab. □for the definition of the energy bands). 


the west is less obvious because of the soft X-ray extinction likely 
caused by a series of molecular clouds (e.g., the Brick) present at 
that location and because of the partial high-latitude coverage of 
this region. 


8.7.2 Radiative process 

The presence of intense, soft X-ray emission lines (see Figs.and 
ED and, in particular, the good fit of a spectral decomposition based 
on a thermally emitting gas (see Fig. ED indicate that most of the 
high-latitude emission is generated by a thermal radiative process 
in a warm plasma. To demonstrate this, we accumulate the EPIC- 
pn spectrum from a circular region of 8.28 arcmin radius centered 
at / = 0.181°, b = 0.359°. The resulting spectrum is well fitted 
with an absorbed thermal emission component (APEC) with kT — 
0.96 ± 0.1 keV, Nh = (2.3 ± 0.2) x 10^^ cm-^ and Aapec = 
(1.5 ±0.4) X 10“2 cm“®. 

8.7.3 Eastern edge 

The eastern edge of the high-latitude emission rises from the po¬ 
sition of the Radio Arc. This raises the interesting question of 
whether the soft X-ray emission at the location of the Radio Arc 


might have two contributions, one associated with the GO. 1-0.1 
superbubble (filling the mid-IR arc bubble; see § |8.6| ), while the 
second is associated with enhanced soft X-ray emission due to the 
presence of the Radio Arc and its polarized radio plumes at higher 
latitudes (Seiradakis et al. 1985; Tsuboi et al. 1986; Yusef-Zadeh 
& Morris 1988). If that is indeed the case and if the two structures 
have different X-ray colours (e.g. the superbubble produces lower 
temperature thermal X-ray lines, while the Radio Arc has a larger 
continuum to lines ratio), then we should observe variations in the 
X-ray colour distribution. In particular, we would expect a whiter 
colour and a green-yellow colour (similar to the one characteris¬ 
ing the lobes of Sgr A) at the location of the Radio Arc region 
compared to GO. 1-0.1, in the top panel of Fig. and in Fig. ED 
respectively. This idea is, indeed, in agreement with the colour vari¬ 
ations and the evolution of the line intensities observed between the 
GO. 1-0.1 and Radio Arc complexes (Fig.p^andpYj). 

8.7.4 Western edge, the Chimney andAFGL5376 

Running almost parallel to the western edge of the high-latitude 
plasma is another region of enhanced soft X-ray emission, located 
near the Galactic plane, the so called Chimney (1 = 359.45°; Tsuru 
et al. 2009). The Chimney appears as a column of soft X-ray emit¬ 
ting plasma extending all the way between the core of Sgr C and 
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Figure 20. (Left panel) 20 pirn MSX map of the GC. The contours indicate the intensity of S xv emission. Soft X-ray emission fills the arc bubble observed in 
the mid-IR. The green solid circle and the white dashed ellipses indicate the position of the PWN GO. 13-0.11 and three structures in the soft X-ray emission 
map (see right panel). (Right panel) Soft X-ray map of the GC (the same energy bands used in Fig.Q^are displayed). The position of the PWN GO.13-0.11 
is indicated by a red circle with 1.5 arcmin radius. At least three sub-structures, appearing like holes, are observed within GO. 1-0.1 (here indicated with white 
dashed ellipses). The positions of the Arches and Quintuplet star clusters are indicated by yellow dashed circles. The direction of the supersonic motion of 
the Quintuplet cluster is indicated and its past location is indicated by the yellow dashed line. The inferred positions of the Quintuplet 4x10^ and 9x10^ 
years ago are indicated with yellow crosses. The cyan dashed ellipse indicates the region in which many massive stars that might have been expelled by the 
Quintuplet cluster are located. 


the northern limit of the XMM-Newton scan (b ~ 0.15°; see Fig. 

pT] and [^. Tsuru et al. (2009) suggested that the Chimney 
is an outflow emanating from the supernova remnant candidate 
G359.41-0.12. They estimated a thermal energy and dynamical 
time for G359.41-0.12 and the Chimney of Eth = 5.9 x 10^® erg, 
tdy = 2.4 X 10^ yr and Eth = 7.6 x 10^^ erg, = 4 x 10^ yr, 
respectively. The energetics and time-scales are consistent with typ¬ 
ical GC supernova remnants and Tsuru et al. suggested that the very 
peculiar morphology of the outflow producing the Chimney might 
be due to a peculiar distribution of molecular clouds that block the 
plasma expansion in the other directions (Tsuru et al. 2009). We 
note that the morphology of the Chimney resembles that of the 
Radio Arc. It originates near the Galactic plane (where dense and 
massive molecular clouds are located) and extends almost perpen¬ 
dicular to the Galactic plane. Within the gap between the Chimney 
and the western edge (see Fig. 6), a bright non-thermal radio fil¬ 
ament with an X-ray counterpart is observed: G359.54-i-0.18 - the 
Ripple filament, with a radio length of 0.08° (Lu et al. 2003; Yusef- 
Zadeh et al. 2005). It is oriented parallel to the edge of the soft X- 
ray plasma distribution (Bally et al. 1989; Yusef-Zadeh et al. 1997; 
2004; Staguhn et al. 1998; see also Yamauchi et al. 2014). Simi¬ 
lar to the Radio Arc, other X-ray and non-thermal radio filaments 
are observed at the base of the Chimney. The high concentration of 
non-thermal filaments indicates the importance of magnetic struc¬ 
tures in this region (e.g., Morris 2014). 

The bright IR source AFGL5376 is located further to the 
northwest, along the continuation of the sharp X-ray edge and the 
Chimney (unfortunately just off ih^ XMM-Newton map; see Uchida 
et al. 1990; 1994). It is associated with high-velocity CO emis¬ 
sion and defines the most prominent portion of a strong large scale 
(~ 90 pc) shock front that extends all the way down to Sgr C 
(Uchida et al. 1994). Because the Chimney appears to be associated 
with a shock, and because it is spatially coincident with magnetic 


filaments along its length (Yusef-Zadeh et al. 2004), we suggest 
that it is not a simple supernova remnant, but is a phenomenon as¬ 
sociated with a footpoint of a larger scale structure, the GCL. 

8.7.5 Is the outflow conflned inside the GCL? Inside the Fermi 
bubbles? 

The GC is considered a mini-starburst environment, producing in¬ 
tense outflows (Crocker 2012; Yoast-Hull et al. 2014). The warm 
plasma detected at high-latitudes is therefore, most probably, as¬ 
sociated with intense star formation and it can be a pervasive at¬ 
mosphere above the entire CMZ. Even in the absence of another 
confining force, the gravitational potential (e.g., that of Breitschw- 
erdt et al. 1991; Launhardt et al. 2002) would bind the ~ 1 keV 
plasma (having a sound speed of ~ 500 km s“^; Muno et al. 2004) 
to the Galaxy, but in hydrostatic equilibrium, would allow it to ex¬ 
tend to heights of several hundred parsecs. If so, it would require a 
large average star formation rate and concomitant energy input to 
generate and maintain it. 

However, the detection of edges in its distribution suggests 
that such plasma might be conflned within known structures. As 
noted by Blanton (2008), the locations of the eastern and western 
footpoints of the GCL (Law 2011) correspond to the positions of 
the Radio Arc and the Chimney, respectively. Indeed, the GCL and 
its possible magnetic nature might confine the warm plasma ob¬ 
served in soft X-rays. This opens the exciting possibility that the 
observed high-latitude enhanced X-ray emission from the GC ’’at¬ 
mosphere” is indeed the warm plasma Ailing the GCL. Based on 
the spectral fit, we deduce a density of rie = 0.06 cm“^ inside 
the GCL. Assuming uniform physical conditions inside the GCL (a 
cylinder of 45 pc radius and 160 pc height) and extrapolating over 
the entire GCL, we estimate a mass of ~ 4 X 10“^ Mq filling the 
GCL with a total thermal energy of ~ 10®^ erg. This value is of the 
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same order of magnitude as the energy required to inflate the GCL 
as estimated by Law (2011). 

Just after the first detection of the GCL, Uchida et al. (1985) 
noted its similarity with the lobes in nearby radio galaxies (al¬ 
though smaller in size and strength). The authors interpreted the 
lobes as created by a magneto-dynamic acceleration mechanism 
where the magnetic twist is produced by the rotation of a contract¬ 
ing disc of gas in the Galactic plane. Under such conditions, the 
plasma is accelerated into a conical cylinder with a helical velocity 
held (Uchida et al. 1985). Alternatively, Bland-Hawthorn & Co¬ 
hen (2003) suggested that the GCL could be produced by a large- 
scale bipolar Galactic wind, that would be the result of a power¬ 
ful {E = ]^q 54-55 nuclear starburst that took place a few 
10® yr ago. These authors show that dust is associated with the 
entire GCL structure and they suggest that the GC (and the cen¬ 
ters of many Galaxies) would drive large-scale winds into the halo 
with a recurrence time of about 10 Myr (Bland-Hawthorn & Cohen 
2003). Other, alternative, scenarios for the origin of the GCL in¬ 
volve outflows associated with enhanced activity of Sgr A* (Ponti 
et al. 2013) or intense star formation (Crocker et al. 2011; 2012). 

It is not excluded that the GCL could be simply one part of 
an even larger scale feature extending over a physical scale of sev¬ 
eral kiloparsecs above and below the Galaxy, the so called Fermi 
bubbles (Su et al. 2010). These gamma-ray bubbles, detected with 
Fermi, are interpreted as produced by highly relativistic particles 
emitting brightly at GeV energies and beyond and they appear 
to contain and confine soft X-ray emitting plasma traced by the 
ROSAT all sky survey, from very large scales down to the Milky 
Way’s center (Su et al. 2010). However, close to the Galactic plane, 
the bubbles’ edges start to become confused. Whatever their origin 
might be, the Fermi bubbles appear to originate (and be collimated) 
from the CMZ, within the region that XMM-Newton scanned here. 

Additional XMM-Newton observations at high Galactic lati¬ 
tudes, in particular, inside and at the border of the GCL, cover¬ 
ing the AFGL 5376 source and the edges of the base of the Fermi 
bubbles will be needed to measure the extent of this high-latitude 
emission and to help disentangle the hypotheses for its origin. Fur¬ 
thermore, higher spatial resolution observations (such as provided 
by Chandra) at high latitudes would allow one to pin down what 
fraction of the extended, high-latitude X-ray emission is associated 
with faint point sources that are relatively less subject to extinction 
than sources near the plane. 

8.8 Soft X-ray emission from the Sgr D and Sgr E regions 

Intense, soft, diffuse X-ray emission is observed from G359.12- 
0.05, the region around 1E1740.7-2942. A radio SNR (G359.07- 
0.02) is observed at about the same position (LaRosa et al. 
2000). G359.12-0.05 has an emission spectrum typical of an SNR 
(Nakashima et al. 2010). In particular, the high extinction sug¬ 
gests it is located at the GC. Nakashima et al. (2010) suggest 
that G359.12-0.05 might be associated with the great annihilator 
and therefore be the second system (such as SS433 and the radio 
SNR W 50) where a BH is associated with its SNR. 

The core of the Sgr D complex is also observed to show en¬ 
hanced medium energy emission (see Fig. [^. In radio, a SNR 
southwest of Sgr D’s core and H II regions are clearly observed 
(Fig.[^ LaRosa et al. 2000). Sawada et al. (2009) analysed X-ray 
data from the XMM-Newton and Suzaku satellites and observed soft 
X-ray emission from two diffuse X-ray sources, DSl and one asso¬ 
ciated with the core of Sgr D. They suggest that DSl is a new SNR 
at the GC. 


8.9 Star formation estimate from counts of supernova 
remnants 

We observe a total of ~ 10 — 12 supernova remnant candidates 
in the CMZ (plus ~ 5 independent radio SNR; see Tab. plus 
two superbubbles, each likely created by many (3-10) supernova 
events. These remnants have typical estimated ages of a few tens 
of thousands of years and temperatures of kT ~ 0.4 — 1.5 keV. 
Due to the presence of the two superbubbles requiring multiple 
SNe and the high absorption towards the GC (such as in the star 
forming region Sgr B) that hampers us from observing a poten¬ 
tially larger population of remnants characterised by lower tem¬ 
peratures, the number of SNR observed in the GC is most prob¬ 
ably under-estimated. However, assuming lifetimes of 10-40 kyr, 
the observed number of SNR yields a rate, averaged over the past 
several thousands of years, of rsN ~ 3.5 — 15 x 10“^ yr~^, 
consistent with other estimates (Crocker 2012). This implies a ki¬ 
netic energy input higher than 1.1 x 10^® erg s“^. To estimate the 
star formation rate, we assume that all stars with masses greater 
than 8 M© produce supernovae and that all SNR are observable. 
Therefore, we multiply the supernova rate by the integral of the 
initial mass function (IMF) over all masses divided by the integral 
of the IMF above 8 M©. To reflect the GC IMF, we assume the 
Kroupa (2002) formulation. The star formation rate then results to 
be: rsFR ~ 0.035 — 0.15 M© yr“^. If the IMF in the CMZ is 
top-heavy, as some have argued, then a smaller star formation rate 
is implied. 

As noted also by Mori et al. (2008; 2009) and Heard & 
Warwick (2013), we observe that the two superbubbles have far 
hotter temperatures (higher density and smaller size) than all the 
ones observed in the Galactic plane or in the Large Magellanic 
Cloud (typically with temperatures kT ~ 0.1 — 0.3 keV, densi¬ 
ties He ~ 0.01 — 0.03 cm“^, sizes I ~ 140 — 450 pc; but see 
also Sasaki et al. 2011; Kavanagh et al. 2012). This could simply 
be the consequence of the high extinction towards the GC, hiding 
a population of normal and very soft superbubbles, or it could be a 
characteristic feature of GC superbubbles, inducing a different evo¬ 
lution because of the interaction with the peculiar GC environment. 
Further investigation is required to solve this problem. 


9 CONCLUSIONS 

We have systematically analysed more than 100 XMM-Newton ob¬ 
servations pointed within one degree of Sgr A* and have created 
the deepest, few arcsec resolution. X-ray images of the CMZ. This 
includes a total of about 1.5 Ms of EPIC-pn cleaned exposure in the 
central 15 arcmin and about 200 ks at all other points of the Cen¬ 
tral Molecular Zone (CMZ). We present here, for the first time, not 
only broad-band X-ray continuum maps, but also mosaicked maps 
of both soft line intensities and inter-line emission from the entire 
CMZ region. 

• The remarkably similar distributions of both the soft line emit¬ 
ting plasma (Si Xlll, S XV, Ar XVII and Ca xix) and the soft contin¬ 
uum (intra-line bands) indicate that most of the diffuse soft X-ray 
emission arises from a thermal process generating both continuum 
and lines. 

• Starting from the mosaic maps of the different narrow en¬ 
ergy bands and assuming the GC emission is produced by three 
different components, we fit the maps at different energies and 
derive the integrated intensity map of the thermal plasma emis¬ 
sion. Integrating over the entire CMZ, the total observed (un- 
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absorbed) flux is: F 2 - 4.5 kev — 4.2 x 10“^^ erg cm“^ s“^, 
Fa .^-12 keV = 1.2 X 10“^° erg cm“^ s“^, corresponding to a 
luminosity of L 2-12 = 3.4 x 10^® erg s“^ at an assumed 8 kpc 
distance. 

• Counting the number of supernova remnants in the CMZ, we 
estimate a supernova rate between rsN ~ 3.5 — 15 x 10“^ yr“^, 
consistent with other estimates (Crocker 2012), that corresponds 
to a star formation rate of tsfr ~ 0.035 — 0.15 M© yr“^ over 
the past several thousand years. This implies a kinetic energy input 
greater than 1.5 x 10^° erg s“^. 

• We report the discovery of a new X-ray filament 
XMM JO.173-0.413 perpendicular to the Galactic plane and 
south of the GC Radio Arc spatially corresponding to a non- 
thermal radio filament. XMM JO. 173-0.413 is the first of the four 
cases known where the X-ray emission is not at or near a location 
where the radio filaments show unusually strong curvature. 

• The soft GC X-ray emission is absorbed not only by high- 
column-density foreground clouds located in the Galactic disk, but 
also by some clouds located on the near side of the CMZ, such as 
the core and envelope of Sgr B2, M0.25-i-0.01 (the ’’Brick”), and 
even a few clouds at higher Galactic latitudes, MO. 18-1-0.126 and 
MO.20-0.48. However, the majority of the observed variations in 
the soft X-ray emission are true emissivity modulations and not a 
product of absorption. 

• Several SNR candidates are identified by their soft X-ray 
emission that appears to All holes in the column density distribu¬ 
tion of gas-and-dust derived from submillimeter maps. 

• Our data shed new light on two quasi-symmetric lobes situated 
to Galactic north and south of Sgr A*. The Northern lobe shows a 
bright and sharp transition at its edge, suggesting the presence of a 
shock. Such features are possibly the remnant of the SN that gen¬ 
erated SGR J1745-2900 or the PWN candidate G359.945-0.044. 
Alternatively, the lobes might constitute a long-lived bipolar struc¬ 
ture produced by an isotropic outflow produced by either 1) the 
cumulative winds from the young stars of the central cluster, 2) a 
wind associated with the accretion flow onto Sgr A*, or 3) the same 
process that generated the X-ray reflection nebulae (if such activity 
has been recurrent over the past millennia). 

• The uniform X-ray colour of the superbubble G359.9-0.125, 
its sharp external edge and its being filled with S XV emit¬ 
ting plasma suggest that the soft X-ray features southwest of 
Sgr A* form a unique shell-like structure with total energy Eth ^ 
10^^ erg, therefore making it a superbubble candidate in the GC 
(high absorption indicates that G359.9-0.125 is located at the GC). 
Alternatively, it might be the remnant of a very energetic event at 
the GC, such as a tidal disruption event. 

• We discover new evidence for the GC superbubble GO. 1-0.1, 
also known as the arc-bubble from mid-IR observations: its soft X- 
ray (e.g. S XV) emission completely Alls the mid-IR bubble, and 
indicates a thermal energy as large as Etn ~ 1-5 x 10^^ erg. At 
present the Quintuplet cluster, which is moving at very high speed 
through the CMZ, is located at the border of the superbubble. How¬ 
ever, it was more centrally located a few 10^ yr ago and it could 
have, at least in part, energised it. We do not observe similar soft X- 
ray emission trailing the Arches cluster, but this might be ascribed 
to its younger age. 

• We suggest that the Galactic Center Lobe might be a magnetic 
structure filled with warm, soft. X-ray-emitting plasma. In fact, we 
observe: i) enhanced soft X-ray emission at high Galactic latitudes; 
ii) enhanced soft X-ray emission at, and between, the longitudes 
of the Radio Arc and the Chimney associated with Sgr C, corre¬ 
sponding to the east and west foot-points of the GCL; iii) a sharp 


edge (at I = 359.63°, 5 = 0.06° and I = 359.55°, 5 = 0.46°), 
running parallel to the nonthermal ripple filament (G359.54-i-0.18) 
and Sgr C thread, defining the western border of the enhanced soft 
X-ray emission. The GCL could be the relatively small base of an 
even larger structure, the so-called Fermi Bubbles. Additional ob¬ 
servations will be needed to clarify this. 

• A new very faint X-ray transient, XMMU J17450.3-291445, 
has been discovered during the new XMM-Newton campaign to 
reach a peak luminosity of Lx ~ 10^^ erg s“^ for ~ 2 hr (Soldi et 
al. 2014). 
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Table 8. List of all XMM-Newton observations considered in this work. Exposure Mode: U, S stand for unscheduled and scheduled, respectively. Filters: Med, 
thn, tck stand for medium, thin and thick filters, respectively. FF, eFF, SW, Ti, TU stand for full frame, extended full frame, small window, timing and time 
uncompressed, respectively. 
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Table 9. List of all XMM-Newton observations considered in this work. Exposure Mode: U, S stand for unscheduled and scheduled, respectively. Filters: Med, 
thn, tck stand for medium, thin and thick filters, respectively. FF, eFF, SW, Ti, TU stand for full frame, extended full frame, small window, timing and time 
uncompressed, respectively. 
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